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Abstract 

We present a refined and expanded analysis of tlie CDF ee^j + flx event as superpartner 
production, assuming the lightest neutralino is the lightest supersymmetric particle. A 
general low-energy Lagrangian is constrained by a minimum cross section times branching 
ratio into two electrons and two photons, kinematics consistent with the event, and LEPl- 
LEP130 data. We examine how the supersymmetric parameters depend on the kinematics, 
branching ratios and experimental predictions with a selectron interpretation of the event, 
and discuss to what extent these are modified by other interpretations. Predictions for 
imminent CERN LEP upgrades and the present and future Fermilab Tevatron are presented. 
Finally, we briefly discuss the possible connection to other phenomena including a light stop, 
the neutralino relic density, the shift in Rh and the associated shift in as, and implications 
for the form of the theory. 
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1 Introduction 



Minimal low energy supersymmetry provides the most promising framework to extend 
the Standard Model (SM). Such extensions take the form of complete models that encompass 
the gauge group structure and particle content of the SM, along with the supersymmetrized 
interactions and superpartners. General low energy theories of supersymmetry have over 100 
parameters in addition to the SM parameters; such parameters can certainly be constrained by 
direct collider searches, but in general one needs more information or more assumptions to do 
calculations that examine many parts of the remaining parameter space. In many cases only 
one or a few parameters enter the calculation of a given observable, so useful predictions can 
often be made from a small subset of the super symmetric parameters without loss of generality. 
The two obvious approaches to reduce the parameter space are to use theoretical assumptions, 
and (direct and indirect) experimental constraints. 

In Ref. [H we showed that the CDF 6677 + If^T event ||2| at the Fermilab Tevatron could be 
interpreted in low energy supersymmetry with roughly the expected rate and kinematics. If we 
assume this interpretation is correct and the event is due to supersymmetry, then we can reduce 
the parameter space by searching for sets of parameters that satisfy the event's constraints. We 
use the term 'model' to describe a distinct set of parameters, but of course all of our 'models' 
parameterize only one basic supersymmetric low energy Lagrangian. The primary difficulty in 
deriving precise parameter constraints (hence predictions) is the somewhat arbitrary notion of 
interpreting one event in terms of a cross section times branching ratio. Instead of advocating 
a particular lower (or upper) threshold value, we vary the value in a reasonable range and show 
the effect on parameter space and predictions. In this way we attempt to give an appreciation 
for the robustness or confidence of particular constraints or predictions. 

We work within a general low energy (= electroweak scale) supersymmetric theory without 
assuming common scalar or gaugino masses at the unification scale |^] . To determine branching 
ratios and scalar interaction contributions to cross sections, we do assume squark mass degener- 
acy except possibly for the light stop ti, and a mass degeneracy among sleptons with the same 
electroweak quantum numbers. Such assumptions are not crucial to our analysis, and could be 
removed if necessary. We assume R-parity is exactly conserved, so the lightest supersymmetric 
particle (LSP) is stable (consistent with the 6677 + fx event where the two LSPs escape the 
CDF detector) . Finally, throughout this paper we assume the LSP is the lightest neutralino A^i , 
and not the gravitino. Analyses of the 6677 + event assuming the LSP is a light gravitino 
have been presented by us ||l|, and in other Refs. [^, |6|. One cannot distinguish these scenarios 
based solely on the 6677 + fx event, although it is likely that associated phenomenology can 
distinguish the scenarios. In this paper we assume that A^i is the LSP, or is at least long-lived 
enough to escape the detector. If Ni is identified as a stable LSP, then it is a possible cold dark 
matter particle Q. 
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In minimal low energy supersymmetry the possibility of one-loop radiative decay of neu- 
tralinos |Q, ^, |l^, |l^ leads to signals with hard isolated photons plus missing energy in the final 
state, a signal predicted many years prior to the 6677 + event. This is by no means the 
only mechanism to produce photons plus missing energy, but it does allow the interpretation 
of the 6677 + $T event as selectron production pp — > e+e~(-|-X), with the selectron e decaying 
mainly into the next-to-lightest neutralino A''2 and an electron, followed by A''2 — > iVi7- It is 
also possible to imagine other interpretations that involve the radiative decay of but for 
which the initial superpartner production is different. The two possibilities in this class that 
we consider below are chargino pair production and neutralino pair production. 

The plan of the paper is as follows. In Sec. |2| we discuss the kinematics of the 6677 -|- IpT 
event in the selectron interpretation, the chargino interpretation, the neutralino interpretation, 
and other interpretations. Using superpartner mass constraints established from the ee^^ + ]pT 
event kinematics, we discuss low energy supersymmetric model building in Sec. |. Here we 
present a discussion of the radiative neutralino branching ratio, slepton decay and constraints 
from LEP. In Sec. ^ we discuss the results obtained from a numerical scan of the parameter 
space, using the structure built up from Sec. |3[ The bulk of our results are contained in Sec. |^, 
where we discuss the model building results, the chargino/neutralino/slepton branching ratios, 
and predictions for LEP and Tevatron. In Sec. ^ we discuss the possibility of explaining the 
ee'-ff + ^T event with the further assumption of a light stop ii. Finally, in Sec. ^, we present our 
concluding remarks, including a summary of such questions as distinguishing left- and right- 
selectrons, and the main channels that can confirm the ee^^+^T event is due to supersymmetry 
with an LSP=-/Vi. In Appendix ^ we discuss the viability of the chargino interpretation, and 
the results of attempts at model building. In Appendix ^ we give four sample models in the 
selectron interpretation. 

Note added: As we were completing this paper, three other papers appeared which discuss 



the CDF 6677 + event in various contexts |31, B2, 33]. 



2 Kinematics of the 6677 + event 



The kinematical requirements on the intermediate particles involved in the 6677 -|- event 
are stringent, and for completeness we present a refined analysis based on the procedure outlined 
in Ref. Q. There are three basic possibilities for intermediate (s)particles; we will present these 
in terms of LSP= A''i interpretations, but the analysis is generic and could be applied to any set 
of intermediate particles that satisfy the criteria below. All decays are assumed to occur close 
to the apparent vertex, which would be true of any LSP= A'^i interpretation. The procedure we 
use to find kinematical constraints is to begin with the information on the observed particles |^ , 
assume two- or three- body decays as appropriate, randomly select unconstrained momentum 
components of the unobserved particles on both sides of the decay chain, and then reconstruct 
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the intermediate particle masses based on all possible pairings of electrons and photons. The 
masses of identical particles on both sides of the decay chain are required to be within 2.5 GeV 
to 'pass' the kinematic cut. The net transverse momentum in the event from adding both the 
observed particles and the LSPs is assumed to be \pt\ ^ 20 GeV. 

2.1 Selectron interpretation 

The first possibility is selectron production pp e~^e~{+X) and decay via the 2-body 
mode e — > eA^2 followed by N2 Ni'j. All sparticles are assumed to be on mass shell. The 
general result is summarized in Fig. |l|, where the allowed regions in the mg-mp^^ plane are 
given for a series of maximum values of rrij:^^. The choice to cut off the graph at rrig = 140 
GeV is motivated by a rough lower limit on the selectron cross section, which will be made 
precise in Sec. [4.2[ Since the electron and photon momenta have experimental uncertainties, 
the kinematic results that we derive from the event will have associated uncertainties. Analytic 
forms of the constraints have been extracted and are presented in Table |l|; a few observations 
are in order that will be useful in model building: 

1. m^^ < (50,74) GeV, for mg < (115,137) GeV. 

2. — > 21 GeV, this value increasing to 30 GeV GeV. 

3. rrie — nifj^ > 20 GeV, this value increasing for decreasing rrig. 

4. Given m^^ ^ 33 GeV, then mg ^ 100 GeV. 

5. Only one pairing of electron and photon gives consistent kinematics for ^ 125 GeV. 

The non-trivial mass differences that are required are not surprising, since all of the parti- 
cles in the event have large (transverse) energy. We incorporate the mass difference constraints 
as well as the constraints on the ranges of m^^, rrij^^, and mg in our model building efforts. 

2.2 Chargino interpretation 

The second possibility is chargino production pp CiCj = 1,2), with three possible 
decay chains: 3-body C N2eVe (through an off-shell or possibly on-shell VF), 2-body C 
eVe or C — > t'eC. For either 2-body decay, the on-shell slepton proceeds through another 2- 
body decay e(t'e) — > eit^e)^2, then the photons are obtained through N2 Nij. Calculating 
consistent kinematics requires specifying the six unknown momenta of the two neutrinos as 
well as the unknown LSP momenta in the final state. This is too complicated to delineate 
any rigorous exclusion regions using the randomized momenta procedure as in the selectron 
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Figure 1: The kinematieally allowed region of the ee^^+J^T event in the mg-mjy^^ plane is shown 
for various values of m^^ in the selectron interpretation. The allowed regions for m^y^ =0, 10, 
20, 30, 40, 50, 60, 70 GeV are to the inside and right of the indicated lines. The allowed region 
for any given m is roughly a subset of any lower m , except for large values of m^^^ . Since 
the lines are derived from the momenta of the 6677 + event, they are only as precise as the 
associated measurement of momenta. 
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Table 1: Kinematical constraints in the selectron interpretation. 
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interpretation. However, we have checked that it is possible to generate consistent kinematics 
for > 95 GeV, assuming the 2-body decay C — > W and that all (s)particles are on-shell. The 
rough regions where we were able to find kinematical solutions have m^- ^ 60, 75 for slepton = 
sneutrino, selectron. In addition, we found solutions with ^ 20 GeV, — ^ 10 
GeV, ^ max[2.5m^^ — 95, 1.5m^^ + 65, 95] GeV. Thus a solution in the selectron 
interpretation need not be a solution in the chargino interpretation, and vice versa. 

2.3 Neutralino interpretation 

The third possibility is neutralino production, e.g. pp N2Nj, where either of the heavier 
neutralinos j = 3,4 decay as Nj l^l^N2, followed by the usual ^i7- This interpretation 

contrasts with the first two by producing both leptons from one side of the decay, however it 
is calculable as in the selectron scenario (since the only unknown final state momenta are the 
two neutralinos). The invariant mass of the electron pair can be extracted from the event 
rrig+g- ~ 160 GeV which implies the mass difference between m^y and must also be 
greater than 160 GeV. This is almost certainly too high for a reasonable Tevatron cross section 
while retaining a reasonable "mj^^ and proper neutralino mixing to have N2 — > Further, in 
the particular case where the branching ratio for the decay Nj — > N2Z is large, then a lepton 
pair from Z will always reconstruct to to an invariant mass of about mz- Thus, a 

neutralino interpretation of the 6677 + event seems extremely unlikely, and we will not 
consider it further. 

2.4 Other interpretations 

Other supersymmetric interpretations with a neutralino LSP are in principle possible, and 
are based on variants of selectron production, chargino production or neutralino production. 
The differences lie in the particular decay from which the electrons originate, plus possibly 
other invisible phenomena (neutrinos). In all cases the photon is obtained from the decay 
N2 -^17) and as a consequence the photon always appears in the last step of the decay 
chain. One example is stau production pp — > f^f~ {+X) with the decay f — > TN2, followed by 
r ^ e (+ i^T-i^e)- The total branching ratio is suppressed compared with selectron production 
by a factor B{t — > eur^e)^ ~ 0.03, hence the rate into 6677 is much smaller than selectron 
production. Another example is a variant of the selectron interpretation with a chargino C 
that is lighter than the selectron, such that the decay ei — > VeC{-^ N2eV(,) is dominant. In 
this case it is probably not possible to have a large decay — > Cve, with both Ni^2e 
suppressed. Further, C — > Nieve has to be suppressed with respect to C — > N2eue, which 
is difficult especially in the presence of N2 — > A^i7. Finally, with four neutrinos carrying off 
invisible momentum it seems difficult to have a large probability for the high energy electrons 
required in the final state, since the selectrons have to be light to have a large 6677 rate. 
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3 Model building 



The kinematics of the event have illustrated two viable sources of ee^^ + ^x events: slepton 
production or chargino production. In either case, the essential ingredient to getting photons 
is through the one-loop radiative decay of neutralinos. To proceed, we first define the relevant 
parameters of the low-energy super symmetric theory, including the chargino and neutralino 
mass matrices. This sets the stage for the discussion of the radiative neutralino branching 
ratio. We also discuss the treatment of the squark, slepton and Higgs sectors and the relevant 
mixings, as well as discussing the selectron branching ratios. Once the models have been 
constructed, we describe the constraints imposed on the parameters from experiment. 

The main focus of this paper is on the selectron interpretation and not the chargino inter- 
pretation, since it is made clear in Appendix ^ that the chargino interpretation is difficult for 
many reasons. However, in the following we have attempted to provide a general discussion of 
the model building, since radiative neutralino decay is required in both interpretations. 

3.1 Supersymmetric parameters 

The chargino and neutralino tree-level masses and mixings are determined by specify- 
ing the gaugino soft masses Mi and M2, the ratio of the Higgs vacuum expectation values 
tan P = {H2)/{Hi) and the Higgs superfield mass parameter /x. The form of the mass matrices 
is well known, but it will prove useful in the discussion of the radiative branching ratio to have 
the expressions in the particular basis as follows. Note that we assume no relation between Mi 
and M2. 

The chargino mass matrix in the {—iW^, H^) basis is 





and can be diagonalized by a biunitary transformation U* M^±V to yield the masses and 
mixing matrices [/, V (as well as fixing the sign convention of consistent with Ref. |12U). The 
chargino masses can be found from the analytic expression 




^i/(M| - iJ?Y + AM^ cos2 2^ + 4M2,(M| + ^2 + 2M2/X sin 2/3) |. (2) 
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Process LEP 



Tevatron 



NiNj meL,meji m^j^, m^^, m^^, mj. 
NjCf - niu^ , m 



Table 2: Chargino and neutralino cross sections at LEP and Tevatron depend on M\, M2, tan /3, 
H and the particular superpartner masses as above. (The Tevatron cross sections also depend on 
the second family masses, but these contributions are generally suppressed by Cabbibo mixing 
and a small parton distribution f^^p in the proton.) 

The neutralino mass matrix in the (— ^7, —iZ, Ha, H^,) basis is 

/ Ml ms^Oyv + M2 sin^Oy^, (M2 - Mi) sin Oy^, cos 6^0 \ 



Mz /i sin 2/3 -;Ucos2/? 

\ -//cos 2/3 -/i sin 2/3 J 



(3) 



and can be diagonalized by a unitary transformation N*Ai ffN^'^ to yield the four neutralino 
mass eigenvalues eirrij^, and the mixing matrix that we assume to be real and orthogonal 
(exact expressions for the mixings and masses can be found in fl^ , p^). The sign of the 
neutralino mass eigenvalue enters the super symmetric Feynman rules, while the physical 
masses are always positive with the ordering < < < ^ri^Vg ^ ^n^- The (7, Z) 
basis is related to the {B, W"^) basis through 




cos sm t>w 
— sin 6t^- cos 6^v 




(4) 



and the {Ha, Hi,) basis is related to the {Hf, H^) basis through 

^ (5) 





Our notation follows Refs. [p!2| , |T^], with H^ and H2 coupling to the down- and up- type fermions 
respectively. The production cross sections for charginos and neutralinos at LEP and at the 
Tevatron involve graphs with s-channel gauge boson exchange and t-channel slepton or squark 
exchange. In Table |2|, we itemize the dependence of each chargino/neutralino cross section on 
the squark or slepton mass. 

The gluino does not enter phenomenology directly associated with the ee^^ + ^T event. Its 
tree-level mass is given by the soft mass parameter M3 that is unconstrained without gaugino 
mass unification. There need be no relation between Mi, M2, and M3, and we do not assume 
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one. However, one could imagine that the non-Abehan masses M2, M3 are equal at the unifica- 
tion scale, with the U(l) mass Mi related to them in a more subtle way. Ref. [^] has suggested 
that the gluino may play a dramatic role at the Tevatron, if the lightest stop ti has a mass 
C(50) GeV. However, for the primary purposes of this paper we can focus on phenomenology 
that is independent of the gluino. In Sec. |^ we elaborate on the possibility of models that can 
generate an 6677 + event with the additional assumption of a light stop. 

The slepton sector is defined by the masses m^-^ and "i^-^, with my related by the 811(2)^ 
sum rule 

ml = mf - M^|cos2/3|, (6) 

for tan/? > 1, and the couplings to gauge bosons and gauginos fixed by the SM gauge group. 
Slepton production cross sections at the Tevatron are given in Refs. [^, ^, and depend only 
on the mass of the slepton. We assume slepton mass degeneracy (motivated by the absence of 
lepton flavor changing decays), although it is not required by the theory nor the 6677 + 
event. Where necessary, we remark on the effect of removing this assumption on associated 
phenomenology. We also assume L-R mixing in the slepton sector can be neglected. 



The squark sector in our model building is defined for simplicity by a common squark mass 
g, the stop masses w-^^, m^^ and the stop mixing angle Oj:. In this way we achieve a useful 
reduction of parameter space through niq = rriu^ = m = rriuj^ = m = • • • , and we further 
assume for simplicity = rriq. These assumptions can be removed if data becomes sensitive 
to them. The stop mass eigenstates are defined by 

/ cos Or sin Or 





(7) 



with the stop trilinear coupling At (and the soft masses rrig, fni^ uniquely determined by ^ 
and the mixing angle 0^, for a given /i and tan (3. We assume all other L-R squark mixing can 
be neglected. 

The Higgs sector is determined from tan /?, the neutral CP-odd Higgs mass m^, and higher 
order corrections 20 1 . We include one-loop corrections from stops [ 2C ] , and neglect all other 



contributions. In this framework we calculate the charged Higgs mass mjj±, the neutral CP- 
even Higgs masses nih, rnn and the mixing angle a from the above parameters. The Higgs sector 
enters the radiative neutralino decay through the charged Higgs boson, and the branching ratios 
for the heavier superpartners into one or more of /i. A, H, or (neglecting off-shell Higgs 
exchange in 3-body C, N ^ C, Nff decays). 



3.2 Radiative decay of neutralinos 



The radiative decay of neutralinos has been well studied |9|, and it suffices to 

review the mechanism that enhances the radiative branching ratio with respect to the traditional 
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3-body ^^"2 Nif f decays, as pertaining to the 6677 + event. We exclusively discuss 
— > -^17) since heavier neutralinos always have sizeable tree- level branching ratios into 2- or 
3-body channels, causing the radiative branching ratio to be negligible. 

There exists both a kinematical and a dynamical mechanism that can give an enhancement 
of the radiative neutralino decay |ll| . The kinematic enhancement can only occur when the 
mass difference m 



N2 



m 



Ni 



is small 0(10) GeV, so that other decay modes are closed or 
suppressed. However, the kinematics in the selectron interpretation enforce ruj^^ — m^^ > 21 
GeV by Observation 2, and so a kinematic enhancement of the radiative branching ratio is not 



crucial for our purposes (although see Sec. 4.2 for exceptions). 



The dynamic enhancement of the radiative decay occurs as follows. First, examine the 
limit when tan/3 — > 1 and (Mi — M2) [^5|; the neutralino mass matrix (already written in 
a suggestive form in Eq. (P)) becomes particularly simple. 



( M2 



V 





M2 

Mz 






Mz 





\ 






for 



tan (3=1 
Ml = M2 



(8) 



In this limit two neutralinos become pure photino (7) and Higgsino (-ff;,) states, with masses 
M2 and respectively. The other two neutralinos are mixtures of Z~Ha, with masses 



M2+^±a/(M2-/x)2+4M| 



(9) 



For pure 7 and Hf, states, the tree- level couplings ^Hi,Z, ^Hbh{A), and Hbff (in the limit rnj ^ 
0) go to zero, leaving the one-loop 'effective' coupling ^Hf,^ dominant. Thus, by associating 
-^1,2 with 7, Hb, then the one loop decay A'^2 -^i7 is dominant. One consequence of requiring 
the two lightest neutralinos to be either of the states 7 or Hb (hence the heavier two neutralino 



masses are given by Eq. (^)) is that the required mass ordering 



m 



Nl,2 



< rrij^^^ implies 



Mi(=M2), |//| < 



M2+^i±J{M2-nf + 4M^ 



(10) 



See Ref. [11] for a more comprehensive treatment of this issue. What is not determined by 
requiring a large radiative branching ratio by this mechanism is which one of the two lightest 
neutralinos is the photino or Higgsino. 

The extent to which a large radiative branching ratio is possible in general (and in particular 
through the dynamical mechanism without the exact relations above) can be evaluated semi- 
analytically and numerically [11|. As an example, Fig. ^(a) shows contours of the branching 
ratio of N2 
GeV, tan /? 



A'^17 in the M1-M2 plane, for ^ 
1.2, and all squarks heavy m^^ = 



-45 GeV, m 



El 



m, 



110 GeV, ruA = 400 
500 GeV. The thick solid line bounding 
the region defined by {NilH^j)'^ {N2\j)'^ > 0.7 anticipates the constraint on selectron decay from 
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Figure 2: (a) Contour plot for the branching ratio of the radiative neutrahno decay 



the Mi~M2 plane for the case tan (3 = 1.2, ^ = —45 GeV, m 



110 GeV, ruf, 



mr 

*1,2 



A^i7 in 
500 



GeV, and tua = 400 GeV. The B{N2 iVn) = 0.9, 0.75, 0.5, 0.3 levels are shown and 
labeled. The LEP excluded region is shaded. The solid thick line outlines the region where 
{Ni\H^)'^ {N2\'y)'^ > 0.7. (b) Contour plot in the same plane with the parameters above, showing 
the mass difference of the two lightest neutralinos in GeV. This figure is a result of the general 



radiative neutrahno decay analysis of Ref. |11] 



the 6677 + f^T event (see Sec. ^]3| below). Contours in the mass difference rrij^j^ — mfj^ > 3, 
10, 20, 40 GeV are shown in Fig. ||(b). Since the selectron interpretation requires a large mass 
difference — > 21 GeV, only a fairly small region of parameter space remains satisfying 
the constraint of a large radiative neutrahno branching ratio. For example, the region bounded 
by B{N2 -/V17) > 0.5, — > 20 GeV, and the LEP exclusion region is characterized by 
roughly 0.6 < M2/M1 < 1.5 for 60 < Mi < 90 GeV, 45 < M2 < 90 GeV, where the constraints 
on M2/M1 are stronger for larger values of Mi, M2. Of course this example only applies to 
the choice of n, tan mg, rriq, niA values as above, but it gives a reasonable illustration of the 
constraints. The region with a large radiative neutrahno decay centered on the line Mi = M2 
persists as is increased or decreased (the region shifts up or down the Mi = M2 line), but 
tends to shrink (and eventually disappear) as tan /? is increased or the squark or slepton masses 
are decreased. 
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3.3 Slepton decay 



In the selection interpretation, the branching ratio of the selectrons e eN2 is crucial to 
produce an 6677 + event. In general, sleptons couple to the gauginos through the usual 
super symmetrized gauge interactions, and also to the Higgsinos through the Yukawa couplings. 
The Yukawa couplings A^- ~ mi /My/ are strongly suppressed by small lepton masses, and for 
our purposes can be neglected. Since the radiative branching ratio ^ iVi7 requires one of 
-^1,2 to be mostly a photino and the other mostly a Higgsino, then the requirement that the 
selectron decays as e — > eN2 implies the photino-Higgsino content of the neutralinos is unique 
and determined 

{Ni\Hbf ^ 1 

(iV2|7)' ~ 1. (11) 
Based on Sec. this implies |;u| < Mi (= M2), in the limit of pure states. 

If the 6677 + $T event is due to eiei production, one must also consider the branching 
fraction of cl to charginos if kinematically accessible. In the kinematics of the selectron inter- 
pretation no such decay was considered, and naively it would seem possible to suppress this 
decay through a judicious choice of chargino mixings. However, it is also possible that ei pro- 
duction occurs with the selectron decay ei CiVe-, then the decay Ci — > eN2- In the tan/3 = 1 
limit (with neutralinos pure states) the masses of the charginos simplifies considerably from 
Eq. (I) to 

' (12) 



1 

Ci,2 2 



M2 + /i ± V(M2 - ^)2 + 



This expression is the same as Eq. (^) with Mz Mw, and shows that the chargino masses 
are directly correlated with the heavier two neutralino masses. It is a simple matter to show 
that m^^^ > ra^^ is always true (in the tan/3 = 1, Mi = M2 limit), while the coupling of Ci to 
-^2 (= 7) and A'^i (= Hi,) is dependent on the gaugino-Higgsino mixings of the chargino. The 
ei^-C-Vf, couplings are also proportional to the gaugino component of C and so a full numerical 
calculation is necessary to determine the relative size of the branching fractions. This will be 



presented in Sec. 4.4 



3.4 Constraints from LEP 



Throughout our analysis, we applied the most updated limits on the supersymmetric pa- 
rameters and bounds on superpartner masses coming from searches at LEPl, as well as the 
more recent run with -y/s = 130.3 and 136.3 GeV (collectively denoted 'LEP130-136') where 
integrated luminosities of about 2.8 and 2.3 pb~^ were accumulated [21|. We also show the 
combined effect of the LEP limits and kinematical constraints on the selectron and light neu- 
tralino masses in the selectron interpretation of the 6677 -|- event, and the derived ranges 
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of /i, Ml and M2 values. The somewhat conservative LEPl bounds we imposed are 



\22., 



^invisible (-^ 


^ SUSY) 


< 


2.3 X 10" 


Artot(^- 


^ SUSY) 


< 


23 MeV 


B{Z - 


N1N2) 


< 


1.2 X 10" 


B{Z - 


N2N2) 


< 


3.5 X 10" 



(13) 



The evaluation of the supersymmetric contribution to the invisible Z width included not only 
the contribution from the direct LSP production Z NiNi, but also the contribution from 
other channels Z Ni{-^ Niui>)Nj{-^ Nivv). These contributions were then subtracted when 
calculating the supersymmetric contributions to the visible Z width. 

The constraints we applied at LEP130-136 are 

C7(e+e" ^ visible SUSY) < 1.8 pb for ^fs = 130.3 GeV 

(T(e+e- ^ visible SUSY) < 2.2 pb for yfs = 136.3 GeV (14) 



corresponding to the 5 visible event level (before detector cuts) for each of the two runs pl| . 
A few remarks on the calculation of the expected total visible supersymmetric cross section 
are in order. First, we considered only the contribution from chargino/neutralino production, 
since charged sleptons relevant to the 6677 + event need to be heavier than 75 GeV just to 
satisfy the kinematics (see Table We require squarks to be heavier than can be produced 
at LEP, except possibly a light stop whose production cross section is always too small to see 
any events at LEP130-136 with the data sample collected. The total visible supersymmetric 
cross section obviously does not include processes like e^e~ — > NiNi, and e~^e~ — > NiNj when 
both Nij Nivi>. This was achieved by doing a complete calculation of the branching ratios 
for chargino/neutralino decays for every model. To ensure the visibility of the signal, we also 
required large enough phase space in the decay of the produced iVj, ^ which in practice 
implied the mass difference rn^j^ — > 10 GeV, in accord with |21]. 

The following observations are useful to understand in some detail how the LEP constraints 
affect our analysis in a general low energy supersymmetric framework (without assuming any 
relation between Mi and M2). Combining the bounds arising from neutralino searches at 
LEP with the need for a next-to- lightest neutralino rrifi^ > 30 GeV from the 6677 + event 
kinematics (see Sec. ^), one finds the "light Higgsino-gaugino window" with Mi, M2, ^ 
and tan/3 « 1 |2^ is excluded. This also implies > 33 GeV, at least for small tan/3. 
Further, given the light Higgsino-gaugino window is excluded for our purposes, only /i < 
survives LEP constraints such that a large radiative neutralino branching ratio is present |1 1| , 
thus we are left with /x < —33 GeV. For tan /3 > 1.3 either the LEP chargino mass bound 
or the direct search for neutralinos begin to exclude regions with small negative irrespective 
of Ml and M2 values. Given a value of /U, one can find rough regions in the M1-M2 plane 
that are allowed by LEP constraints, generally independent of tan (3. In our framework, the 
constraints we listed above exclude Mi < 30 GeV and, for instance, when /x = —45 GeV then 
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Ml < 55 GeV is not allowed if M2 < 20 GeV. The region in M1-M2 space excluded by LEP 
is indicated in Fig. Q for // = —45 GeV, etc. Notice that since the 6677 + event requires a 
suitable slepton decay, then the neutralino contents in Eq. (^) can exclude a comparable region 
(see Sec. O, and in particular Fig. §). In contrast, the requirement (?7ijv2 ""^TVi) ^ GeV of 



Observation 2 in Sec. y combined with the LEP constraints effectively sets a minimum suitable 
value of Ml around 52 GeV for any values of the other parameters. Only weaker bounds on 
M2 can be identified in a similar way. 

In addition to the constraints from chargino and neutralino production, we also imposed 

f 44 

i^h > i ^r, . . 2fo ^ GeV, (15) 
1 58.4 sm (p — a) 

on our models from LEP constraints. Since the inputs to our model building to calculate the 
Higgs sector include niA and tan (3, the above mass bounds impose a constraint on niA and 
higher order corrections from the stop sector. This will be important for the discussion about 
models with a light stop in Sec. |5[ Small tan /? also suffers from possible non-perturbativity 



constraints, that have been discussed recently in e.g. Ref. |Q for the light Higgsino-gaugino 
window that requires small tan/3. However, the constraint is relatively weak (tan/3 > 1.2), 
since as we shall see the allowed region of tan/3 extends up to tan/3 ~ 2.0 2.8. 



4 Numerical results — selectron interpretation 



To ensure a large branching ratio for the decay N2 — > iVi7, pure photino and Higgsino 
lightest neutralinos are sufficient, but not necessary conditions. The extent of the allowable 
impurity determines the character of the models, but that is by no means the only degree of 
freedom. As we have seen, the branching ratios of the sleptons are also determined by the 
gaugino-Higgsino content of the neutralinos and charginos. Further, the allowed sets of masses 
must satisfy the 6677 + event kinematics, and proper experimental constraints are not 
trivial mass exclusions, etc. What we present here are complete low energy models constructed 



using the framework built up in Sec. y using a randomized parameter selection scheme |24], and 
imposing all of the above constraints. 

4.1 Preliminaries 



Interpreting one event as a cross section is a tenuous procedure, although some general 
methodology can be applied. First, we establish a minimum threshold in the Tevatron selectron 
cross section times branching ratio into two electrons and two photons, 

ax 13^ = a{pp e+e") x fe(e ^ N2e)B{N2 ^ ^^17)] ^ > A, (16) 
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where ^ = (cr x ;B^)|inin is the minimum threshold value. Since the choice of the threshold A is 
somewhat arbitrary, we show the effect of increasing the threshold from 5 to 7.5 to 10 fb to give 
at least some indication as to how sensitive the constraints are to the value. Imposing ^ = 20 fb 
excludes all of our models, so there is a non-trivial importance of the precise numerical value 
of the threshold for phenomenology. 

The quantity a x used in the general analysis does not include detector cuts, but we 



have simulated particular models to get indicative efficiencies (see Sec. 4.8). For a detection 
efficiency of 0.2, the lowest threshold cut ^ = 5 fb corresponds to assuming a cut on the effective 
6677 ^^^^ oi s = a X X EFF = 1 fb or 1/10 of an event. Given an expected number of events 
s, the probability of observing exactly n events is from Poisson statistics 

P=-^- (17) 
n! 

For s = 0.1 corresponding to 1 fb cross section at the Tevatron, one still has a 9% chance of 
seeing exactly one event. 

The results are presented assuming a branching ratio into only one family, although it is 
straightforward to compute the total two lepton plus two photon rate including smuon and/or 
stau production. The effect is of course to increase our calculated rate by a factor of 2 or 3. 
(Our results remain unchanged if the threshold A is increased by the same factor.) Note that 
including more than one family is of course crucially dependent on the assumption of slepton 
mass degeneracy. 

In the selectron interpretation there is no a priori requirement of having ei or e/j pro- 
duction. We consider three cases: A selectron interpretation from ei production, where the 
kinematics of the ee'j'j + ^T event must be satisfied for mg^, but must not be satisfied for mg^. 
In this way, crcr production can still give an 6677 signal but the kinematics are not consistent 
with the 6677 + event; hence only the rate from 61,6^ production ought to be considered. 
Second, the opposite scenario with e/j production where the kinematics must be satisfied for 
mgjj but not for mg^. Finally, we consider a set of models with the simultaneous clcl and 
eneji production (denoted 'e^ -|- cr models'), where the kinematics are satisfied for both rngj^ 
and mgj^. The threshold A is applied as follows, 

aL X B\ > A for cl models 

urx B\ > A for Br models (18) 
ai X Bi + aR X Br > A for ei + e.R models, 

where aL,R = (t{PP e:L,R^L,R) ^'^'^ ^l,r = B{eL,R N2e)B{N2 A^i7)- The case of + cr 
models assumes that the contributions to the 6677 cross section from ei and eR_ production 
can be summed, hence the requirement that the kinematics of the event is satisfied for both 
contributions. Further, for cl+sr models we enforce aL,R > 1 fb to avoid the difficulty of one of 
o'L,R X tS'l R being arbitrarily close, but below the threshold A while the other contribution can 
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be very small. In such a case the model could still pass the cut on the sum x B\+aji x B\ > A, 
but would be on the borderline of classification as either an ei, e_R, or cl + en model. We will 
show that this loose requirement on the cross section does not affect our results. Finally, note 
that since B{N2 -/V17) depends in general on both selectron masses mg^ and mg^, then cl, 
Br and cl + cr models can each be considered a distinct class of models. 

We impose no restriction on the squared branching ratio (unlike Ref. nor any re- 
striction on associated phenomenology. In practice, the cut on o" x does provide an effective 
lower limit on the branching ratio based on the largest allowed cross section a, obtained from 
the smallest selectron mass allowed from 6677 + event kinematics. This avoids generat- 
ing a disproportionate number of non-ee77 events from e^ex, production in models, and 
srsr production in sr models. However, we do not constrain possible non-standard visible 
phenomenology from the other selectron. The absence of knowledge of both the experimental 
data and the efficiency of detection of such phenomenology prevents explicitly restricting our 
models in this regard. As an example, slepton mass degeneracy implies the rate for two smuons 
or staus plus two photons is at the same rate as selectrons. But, without a fully analyzed, 
statistically significant sample of two lepton plus two photon events, one cannot use the lack of 
reported events to exclude such a scenario. 

4.2 Model building results 

In Table ^, we present the parameters that enter our analysis common to all selectron 
interpretations, and the relevant ranges. For the ei and cr interpretations, the allowed range 
of rrie is determined by the lower bound from kinematics nig ^ 100 GeV using Observation 4 
in Sec. § (indeed |//| ~ > 33 GeV, from Sec. 0). The upper bound is obtained from the 
minimum threshold in the cross section times branching ratio A. For ^ = 5, 7.5, 10 fb, the 
upper bound on the slepton mass is mg^ < 137, 125, 118 GeV, and mg^j < 115, 105, 97 GeV, 
in the and cr interpretations. Notice that cr models always fail the highest threshold, since 
the cross section never exceeds 10 fb in the allowed mass range. The mass of the other slepton 
that is not the source of the 6677 -|- event (hence 6677 -|- event kinematics do not apply) 
is allowed to take on a much wider mass range 60-500 GeV. For the ei + gr interpretation, 
both sleptons still must be greater than 100 GeV by 6677 -|- event kinematics, but the upper 
limits are somewhat relaxed since each individual rate ctl x B\ or ur x need not be larger 
than the threshold; only the sum must satisfy the a x B^ constraint. 

We have explicitly constructed roughly 2500 models in total, with somewhat more 
models than cr or CL + eR. The results are shown in a series of scatter plots and bar graphs that 
are intended to give the general character of the models. Figures ^, ^ show the distribution 
of all the allowed models in the M1-M2 plane, with groupings of models split up into three 
plots. All of the models pass the 6677 -|- event kinematic cuts for one or both sleptons 
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Parameter 



Range 



Ml, M2, /U, tan/3 randomized throughout allowed range 
rUg = mi^ 250, 500, 1000 GeV 

nii^ > 150 GeV, < rriq 

O-t [-vr,vr] 
ruA 50, 100, 200, 400 GeV 



Table 3: Parameter ranges common to all selectron interpretations with a heavier stop. Models 
with a light stop are discussed in Sec. ^. 

(defined by the model type), and all models pass the minimum threshold cut „4 = 5 fb. In 
Fig. ^, the models are grouped by the type e^, cr, or ei + sr according to which slepton(s) 
passed the 6677 + event kinematic cuts. In Fig. the models are grouped by the rate, 
5 < cj X < 7.5, 7.5 < a X B'^ < 10, and cr x > 10 fb. In Fig. || the models are grouped 
by tan/3 into the (arbitrary) ranges 1 < tan/3 < 1.5, 1.5 < tan/3 < 2, and tan/3 > 2. There are 
perhaps four regions with distinct character, and we will discuss each of them in the following. 

Region 1 defined by roughly 0.8 < M2/M1 < 1.2 represents the anticipated Mi ~ M2 
region. All three types of models e^, and + eR fall into this range, with cr models almost 
contained within the M2/M1 limits. This is the region where the dynamical enhancement of the 
radiative neutralino branching ratio is present, with the limiting case {Mi — M2) — > 0, tan /3 — > 1 
giving the largest value. Hence, the highest a x B'^ can be found in this region, but the rate 
need not be high since the slepton cross section can be low independent of the branching ratio. 
For example, cr models always have ctr x Bj^ <, 8.2 fb with Bj^ ^ 98%, whereas models 
have ai X B\ < 16.2 fb with B^ < 88%. Since the decay — > Cii^e is always present, the 
maximum branching ratio Bf^ is always less than the maximum for Bj^. 

Region 2 defined by roughly M2/M1 > 1.2 is populated with mostly bl models, extending 
barely up to the M2 = 2Mi line near Mi ~ 60 GeV. The reason for the much larger range in 
M2 values for models is a direct consequence of the higher cross section cjl ~ 2.2a r for a 
given slepton mass. With a higher cross section, the total squared branching ratio can be lower, 
which translates into looser restrictions on the radiative neutralino branching ratio. For and 
cr models, the minimum B^ is 25% and 56%, which corresponds to a minimum B{N2 — > Nij) 
of 50% and 75% respectively (when B{e — > N2e) = 100%). Fig. ^ already showed (for a specific 
set of /i, tan/3, rriq, ruA values) that a looser restriction on the radiative branching ratio admits 
a larger region in the M1-M2 plane. The models observed with M2/M1 > 1.2 lie in just this 
extended region which benefit from the kinematical mechanism (in addition to the dynamical 
mechanism) for the radiative neutralino decay enhancement. This can be deduced by examining 
the slepton masses for the cl models in this region, where one finds m^^ 3> mg^ by a factor of 
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40 I ^ < ^ ^ < ^ < ^ < ^ 1 

50 60 70 80 90 100 110 

Mi [GeV] 



Figure 3: The models satisfying the 6677 + event kinematics and the minimum threshold 
cut ^ = 5 fb are shown in the M1-M2 plane. In this figure, (L), (R) and ei + eji (L+R) 
have been separated to show the varying restrictions on either type of model. 



■ 5<oxB''<7.5fb 




40 I . , ^ . < . < . < . 1 
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Figure 4: As in Fig. ^, except the models are distinguished by the cut on ^ = 5, 7.5, 10 fb. 



18 



140 
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Figure 5: As in Fig. ^ except that the models are distinguished by the value of tan p. 




-60 -50 
1^ [GeV] 



Figure 6: All of the models are shown in the /i~tan/? plane, distinguished by the cut on ^ = 5, 
7.5, 10 fb. 
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2 or more. This is necessary to obtain a large radiative neutralino decay, since the branching 
ratio for 3-body decays — ^ Nil^l~ through sleptons cannot be reduced to zero when the 
kinematical mechanism for a large B{N2 -/V17) operates [^]. In addition, the squark masses 
must also be heavy to prevent the analogous 3-body decays mediated by squarks, although the 
choice of rriq > 250 GeV in our models is sufficient. Finally, the existence of only cl models in 
this region is due to the fact that kinematical enhancement of the radiative neutralino decay 
cannot be maximized simultaneously with the — > 21 GeV, and so B"^ cannot be 
very large. Thus, one needs a large cross section to supplement a lower B^, which can only be 
achieved with ei models. 

The character of the 'extended' cl models in Region 2 is more clearly visible in Fig. where 
all of the models have been plotted in the M1-M2 plane distinguished only by the tan [5 value. 
The models with M2/M1 ^ 1.2 always have 1.5 ^ tan/3 ^ 2.8, where the upper limit in 
tan/3 (and Afi, M2) is established by the smallest allowed radiative neutralino branching ratio. 
Indeed, the kinematical mechanism that contributes to the enhanced radiative neutralino decay 
in this region does not necessarily require tan /? ~ 1 |]ll| . In Fig. |^ it is clear that increasing 
the threshold A to 7.5, 10 pb restricts M2/M1 < 1.9, 1.2, and so the existence of models with 
M2 = 2Mi is sensitive to the choice of the minimum threshold. Further, while M2 = 2Mi 
seemingly admits gaugino mass unification, we noted above that for the extended cl models 
rrigj^ ^ mgj^ . Hence, scalar mass unification probably cannot be achieved, at least in the slepton 
sector, and a completely unified scenario seems not to be compatible with the ee-^^ + ^T event. 

In Region 3 loosely defined as M2/M1 0.8, cr models appear near Mi ~ 75 GeV and 
M2 ~ 50 GeV. These models have aR x 5.5 fb and tan /? ~ 2. This is the only region 

where the usual mass hierarchy < M2 can be slightly violated. On closer inspection one 
finds the chargino mass is about ~ 68 GeV. We found no ei models in this region, due to 
the light chargino that induces a large branching ratio for ei — > Cfe over — > N2e. Also, 
the width for the 3-body decay Nie^e~ turns out to be considerably enhanced when 

the cl is light. Hence the radiative neutralino decay is strongly suppressed in such a case, and 
thus ei models cannot be constructed in Region 3. As tan/3 is increased, the chargino mass 
becomes smaller and thus is excluded by LEP130-136 constraints. Lowering tan/3 decreases the 
radiative neutralino branching ratio, and so is excluded by the a xB^ cut. This localized region 
is basically due to a hybrid of the dynamical and kinematical enhancement of the radiative 
neutralino decay. One can use an argument analogous to that used for Region 2, to observe 
that rrigj^ ^ mg^ in all of the models. The absence of light e/, is a consequence of the kinematical 
mechanism at least partly at work. Thus, these models sit at the edge of exclusion, between a 
multitude of constraints. 

Finally, the voids with no models found for Mi > 85 GeV with M2/M1 < 0.8 or 
M2/M1 > 1.2 are excluded by a low radiative neutralino branching ratio. This behavior 
can be discerned from Fig. |2[ but of course the numerical result here encompasses a full range 
of fi and tan /? values. 
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Naively one might think that + gr models can always be constructed from ei or cr 
models, by simply shifting the other slepton mass such that mg^ ~ '"^e^i- This construction 
always satisfies the 6677 + event kinematics, which are of course invariant under L R. 
Indeed, such a construction can work in the region with a dominant dynamical enhancement of 
the radiative neutralino decay. However, the construction need not work in the region where a 
kinematic enhancement of the radiative neutralino decay occurs, such as in Region 2 populated 
by models. As discussed above, mg^ ^ mg^ in this region which prevented 3-body decays 
N2 Nil~^l~ mediated by Ir to overwhelm the radiative decay N2 -^i7- 

In general, ei + e.R models tend to be constrained similar to cr models, but looser bounds 
on M2/M1 are present and larger Mi values accessible. The region with ei + cr models that 
is devoid of or models, defined as Region 4, has the properties that the a x < 7.5 fb 
and tan/3 < 1.5, while simultaneously cjl x ;B| < 5 fb and o-/? x < 5 fb. For Mi(f« M2) ^ 
90 GeV, larger chargino and neutralino masses are allowed than in either ei or cr models. 
In particular, rrifj^ is near the upper bound from 6677 + event kinematics, so presumably 
values of Mi higher than obtained in ei + cr models are not accessible. As for the size of the 
ee77 rate, the maximum (summed) a x B'^ < 19 fb, so it would appear one does not gain 
more than a factor of about 1.2 over the maximum 6677 rate for ei models alone. Further, 
since cl + cr models enlarge the allowed region of parameter space by reducing the minimum 
CL,i? ^ ^\ Ri one can use the results as an indication of the region resulting from relaxing the 
^ = 5 fb cut in cl or cr models separately. It is clear that cl + e.R models have a distinct 
character separate from ei or eR models. 

In Fig. |6| we show the models in the /i-tan fi plane to completely specify the parameters. 
Three features are worthy of explanation: First, the upper and lower limits on are approxi- 
mately the upper and lower limits on ^jy^, since iVi ~ H},. From Observation 1 in Sec. we 
know the upper limit on is 50, 74 GeV for cr and ei models, and this can be translated 
into rough upper limits on The lower limit on ~ > 33 GeV, and the region devoid 
of models in the upper right-hand corner (larger tan/3, smaller come from a confluence of 



LEPl, LEP130-136 and ee'j'j + constraints as explained in in Section 3.4. For example, the 
LEP constraints on chargino and neutralino production forbid models with |/i| < 40 (50) GeV 
for tan/3 > 1.5 (2), once very small \^\ are excluded by 6677 + event kinematics. 



The final allowed ranges of Mi, M2, fi and the ranges of masses w.^^, "ijVa' "^ATg' "^iV4' 
m^_^, m^^ derived from them are presented in Fig. ^. The effect of imposing a stricter cut 
„4 = 5, 7.5, 10 fb is shown, in addition to the ranges for cr models only. The latter is to give 
an idea of the stronger constraints that exist when a specific origin of the 6677 -|- event is 
assumed. Correlations between a selection of chargino/neutralino masses can be discerned from 
Fig. ^. Sleptons can also have correlations with chargino/neutralino masses, which are relevant 
for the branching ratios. We present these mass ranges in Table For example, notice that 
the mass of the slepton satisfying the 6677 -|- IpT event kinematics always obeys mg > rn^j^ . 

Squarks do not play a large role in our analysis, since they are assumed to be heavier 
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Table 4: Ranges of selected mass differences between cl, Ve-, and eR and ctiargino/neutralinos 
in ei, and eR models. 
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Figure 7: The allowed mass spectrum is shown for all models (shaded bands on the left) and 
for models only (thick solid outline on the right). The increasingly darker shades in the 
left-hand column correspond to increasing stricter cuts on ^ = 5, 7.5, 10 fb. As for tan /3, the 
allowed range in all models is 1.0 < tan/3 < (2.8, 2.6, 1.8) for ^ = 5, 7.5, 10 fb respectively. 
The allowed range of tan /? in e/j models only is 1.0 < tan (3 < 2.0. 

than charginos and neutralinos. However, two effects for a given value of the squark mass 
persist: First, in 3-body decays of neutralinos, the t-channel exchange of squarks can lower the 
branching ratio of ^^^2 — ^ -^17? hence the rate a x B^. Second, the stops enter in the loops of the 
one-loop radiative neutralino decay width (since the Yukawa coupling of Hb to t is significant), 
and also tend to slightly decrease the radiative neutralino branching ratio for lighter m^^ ^ 101 • 
With nig = = 250 GeV, we found no models satisfying the ^ = 5 fb cut, and ei or 
+ (iR models always have a x B'^ < 8 fb. 

The effect of different neutral CP-odd Higgs masses tua is primarily confined to the neu- 
tralino branching ratios, although does enter the one- loop radiative neutralino decay width. 
We find that varying rriA from 50 to 400 GeV does not significantly change the size of the ra- 
diative neutralino branching ratio, hence the a x B^ for the 6677 -|- event. 

4.3 Neutralino composition and branching ratios 

In Fig. ^ we show the (maximum) allowed range of the neutralino composition (iVj|x)^ of 
all of the models in the X = 7; ^1 Hb basis. For a given threshold in a x B^ applied to all 
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Figure 8: As in Fig. ^, except that mass differences between certain charginos and neutralinos 
are shown. 
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Figure 9: The allowed range of all four neutralinos' composition {N\x)'^ in terms of the inter- 
action eigenstates X = 7? Ha, Hf^ is shown for all of the models. The thick solid outline 
corresponds to br models only. Bars that touch the x-axis correspond to a neutralino com- 
ponent that can be lower than 10~^; the absence of a bar for the 7 component of N-^ implies 
(iV3|7)2 < 10^2 foj. all models. 
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models, the minimum radiative neutralino branching ratio is always larger for e/j than for cl 
models. A larger minimum radiative neutralino branching ratio implies the constraints on the 
neutralino composition must be similarly stronger, hence the differing notation for all models 
and models in the plot. We make three observations: First, we find that 

iVi-Ffc, iV2~7, (19) 

so the lightest two neutralinos are composed of exactly the content expected from Eq. (|ll|). To 
a lesser extent, 

iVg-i^a, Ni^Z, (20) 

the heavier two neutralinos turn out have a fairly specific composition as well. This will be 
relevant to the branching ratios and cross sections for associated phenomenology. Second, A^i 
tends to have a much larger Z component than A'2. Third, the required purity of the lightest 
neutralinos in models is significant compared with models, and this is perhaps most easily 
observed by looking at for example the photino content of iVi and ^^^2 in Fig. ^. 

In the following discussion of the branching ratios (and the discussion in subsequent sec- 
tions), we discuss only the distinctions between ei and en models, since the branching ratios 
in e-L + en models are a relatively simple extension of and cr separately. The range of 
branching ratios of N2 are shown in Table ^. In the pure state limit N2 = 7, only the radiative 
channel is open for ^^^2- However, the impurity of N2 (see Fig. |9|) causes other modes to have 
non-negligible branching fractions {N2 is somewhat of a special case since the radiative decay 
branching ratio is required to be large). The possible decays for A^2 in our models are: 
iVi"Z", Ci'W", vv, hi, IrI. We use "Z" and 'W to mean the 3-body decay mediated by 
an on- or off-shell Z and W , plus off-shell sleptons and squarks. The rate for the final states 
"Z"^ l^l^ ,vu^qq and 'W— > lv,qq' are determined roughly by the corresponding SM gauge 
boson branching ratios. The only significant deviation from the SM gauge boson branching 
fractions is modes that involve sleptons, since the ee^^ + event requires at least one slepton 
is light. The presence of some modes depends on the particular class of models; for example 
in e-L models, the mode A'^2 IrI is open if m^-^ < rrifj^. This never happens in eR_ models 
since cr N2e is required to obtain the 6677 + event! The 2-body mode A'^2 is 
open if m/j < nij^^, which happens in cr models and could potentially happen in models. 
However, for ei models one never finds decays N2 vv because the mass splitting between v 
and li^ is never more than about 25 GeV (see Table ||). Since there always must be a large mass 
difference between mg^ and vn^^ from 6677 + If^T event kinematics, then the 2-body mode into 
a sneutrino is always closed. 

The A^3 and branching ratio pattern is progressively more complicated than for the 
lighter neutralino due to possible 2-body decays into sleptons and Higgs bosons. For A^s, there 
are several distinct classes of final states: Ni^2^^Z" , Ci'W, Nih{A), lil, IrI, vv; all other 
possible channels are strongly suppressed or forbidden. For example, for the heavier chargino 
one has > rrifr in all of our models, hence A3 decay into C2 is forbidden. 
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The upper limits on the mass differences rrij^^ — mji;^^ < 60 GeV and mj^^— m^^ < 35 GeV 
in our models are crucial to determining the allowed decays of N3. In particular, the decay 
Nih or NiA will only occur when rrih or rriA < 60 GeV, with constraints from 

LEP that exclude rrih < 44 GeV and the coupling sin^(/3 — a) ^ ^ GeV. The restriction on 
the mass of A from LEP that excludes tha < 22 GeV is considerably weaker than the one on 
nifi, and so decays NiA are always possible with an appropriate choice of rriA (provided 

this does not imply an excluded ruh value). The situation is actually considerably more subtle. 
We find decays into iVi "Z" are not suppressed even if decays into the light Higgs h are open, 
with a maximum B{N-^ Nih) ~ 35% while — > A^i^ decay is closed. However, with low 
tan f3 the mass splitting between rriA and nifi tends to be small for niA ^ 50 GeV, and because 
of the couplings, decays into A typically dominate over h if kinematically accessible. In ei (en) 
models, the decay lil (N^ IrI) is always kinematically forbidden. Thus, it is only 

when the other slepton (//j in ei models, II in models) has a mass rrij < rn^j^ that 2-body 
decays into sleptons can dominate. When kinematically accessible, the branching ratio for the 
2-body decay vv can be ~ 100%, and is always larger than decays into liJ by a factor 

of at least 10. This is due to the larger Z impurity in iVs, i.e. (iVslZ)^ ^ (iV3|7)^, and Eq. Q 
requiring my < mj^ . The 3-body decays into the lightest chargino A''3 Ci 'W depend on 
the chargino mixings, but are always smaller than the 3-body decays N-^ Ni "Z" mainly 
due to phase space. The presence of decays into Ci 'W can suppress the branching ratio for 
decays into A^i"Z" by at most a factor of 2, but even then the branching ratios for are still 
larger into iVi "Z" . Also, decays into N2 are strongly suppressed, because of the particular 
neutralino composition in our models. 

The branching ratios of A'4 are quite intricate, however a few features can be discerned. 
The main possible decays include: Ni^2,3^^Z" , Ci^2'W', ^ih{A), and possibly open 2-body 
modes IlI, IrI, i^i^- Since the mass difference — > 67 GeV and can be as large 

as 100 (120) GeV in cr (ei) models, then the decay A^4 Nih is a prominent possibility if 
kinematically allowed. We find that even if A^4 — > NiA is also open, it is always suppressed to 
of order ~ 5% compared with a much larger Nih mode. This is because the N4 composition 
is roughly inverted with respect to the one, which feeds into the A''4 couplings to the Higgs 
sector. In eji models, the 2-body slepton decay iV4 — > Ijil is always open, and can be ~ 100%. 
In bl models, the decay A'4 — > is typically open, but sometimes can be kinematically 
inaccessible. Note that if both A^4 uv and A'4 — > l^l are accessible, then A'4 — > vv always 
overwhelms A'4 — > l^l by at least a factor of 5 due to the large Z component of A'4 (see 
Fig. |9|) and phase space. Similarly, if none of the 2-body modes are open, then the neutralino 
composition of implies Ni "Z" dominates over all other 3-body decays. 

4.4 Chargino composition and branching ratios 

The chargino composition is determined by the mixing matrices U and V ^ as defined in 
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Final State Kinematic Condition Range in models Range in eji models 



iVi"Z" - ^47 ^26 

Nij - 53 ^ 100 74 ^ 100 

Ci'W" m^2>m^^ ^10 ^2 

IJ + IlI mfj^>mi^ 

IrI + IrI m^^ > rn,;-^ ^3 



Table 5: Ranges of selected A''2 branching ratios (in %) in our models. The notation '— > X' 
denotes a range from less than 1% up to X%. The kinematic condition must be satisfied for 
the mode to be open; no kinematic condition implies the mode always open. Note that the 
2-body decays into sleptons sums over all three families, because of the assumption of slepton 
mass degeneracy. The final state into Nie'^e~~ can be enhanced over that expected from Ni^^Z" 
because of light slepton exchange. 



Final State Kinematic Condition Range in cl models Range in en models 



iVi"Z" - ^99 ^99 

Nih m^^ - m^^ > ^29 ^31 

NiA nij^^^ - rrifj^ > ruA ^66 ^71 

Ci'W" - ^34 ^29 

iVa"^" - ^1.5 ^ 1.5 

IiJ + IlI TTifj^ymi^ - -^22 

IrI + IrI mfj^>mi^ ^99 

w + uv rrifj > nil) — > 99 ^99 



Table 6: Ranges of selected A^3 branching ratios (in %), as in Table |[ 
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Figure 10: Scatter plot of the chargino mixing matrix elements Uu, Vn for all models. The 
narrow band of points indicates the presence of strong constraints in our models from the 
6677 + event. 



Sec. |3.1| . U and V (real and orthogonal in our conventions) can be expressed in terms of two 
independent rotation angles (/>^ (see e.g. Ref. however the Dirac nature of the chargino 

spinors does not allow an intuitive identification of their Wino and Higgsino components. Nev- 
ertheless, in Fig. 10 we present the elements |Viip = cos'^ 0+ vs. = cos^ (/>_ to give a 
sense of the constraints that the 6677 + event imposes on the chargino composition. In the 
limit tan /? — > 1, the chargino mass matrix is symmetric which implies (f)- = and so U = V. 



In Fig. 10, this is the diagonal line where = l^iiP, and note that along this line our 

models lie in the region 0.15 < ~ 0.25, due to the mass hierarchy M2 > Here, one 

can identify Ci as mostly a charged Higgsino. For larger tan/3 values, tends to increase, 

while ^ 0.25 throughout our models. 

The branching ratios of Ci are displayed in Table |^, that assumes m^^ is heavier than both 
charginos as in the discussion below. There are only a few possible channels: A'^i,2'W\ Ili^, 
Further, the 3-body decays into N2^W'^ are always < 5% due to the photino nature of the 
Higgsino nature of Ci and phase space. Thus, the 3-body decays into iVi 'W are the typical 
decay pattern. In cl models /l, t' are always heavier than Ci, thus it is only in e/j models that 
2-body channels into lii^ and Dl can possibly be open. When both are allowed, these 2-body 
decays can sum to a branching ratio of ~ 100% (when summed over three families). 

The branching ratios of C2 are displayed in Table ^. The possible decays include: Ni^2,3 'W , 
Ciff (/ = I, 1^, q), II^, vl, and the Higgs channels Ni^2H^, Cih{A). When only 3-body de- 
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Final State Kinematic Condition Range in models Range in e/j models 



Ni 'W 




95 




100 


^ 100 


N2 ' w 








5 


^ 5 


91 










^ 100 












^ 50 



Table 7: Ranges of selected Ci branching ratios (in %), as in Table ^. m^^ > is assumed 
here. 



Final State Kinematic Condition Range in cl models Range in e/j models 



iVi'W" - ^92 ^ 100 

iVs'W" - ^23 ^17 

iVs'W" - ^0.7 ^ 0.3 

Ci// - ^4 ^1 

^ "^'^ — > 95 — > 69 

Zlz^ mf=, > mr -^52 — > 59 



Table 8: Ranges of selected C2 branching ratios (in %) assuming uia ^ 100 GeV, as in Table |. 
m^^ > is assumed here. 

cays are open, iVi'W dominates over all other decays. However, is roughly 1-5%, and 

can be larger than the decays into N2qq' due to the possible enhancement from light slepton 
exchange in the 3-body decay. The 2-body decay C2 vl summed over three families can 
have a branching ratio up to 95%, when it is the only slepton mode open (the remainder is dis- 
tributed to the 3-body decays as above). When both C2 vl and 6*2 — > /lz^ are simultaneously 
open, the sum can be nearly 100%. Finally, the 2-body decay into NiH^ is also possible when 
"i-_ff± ^ 90 (120) GeV, for eji (cl) models, which requires niA ^ 50 (100) GeV. In addition, 
decays into neutral Higgs bosons are possible when — < mh,A- 



4.5 Sneutrino branching ratios 



In the selectron interpretation, sneutrinos do not directly enter the branching ratios relevant 
for the 6677 + event, however the mass of the sneutrino Ve is necessarily smaller than mg^ 
due to the sum rule in Eq. (P), and so the sneutrino is certainly relevant in and cl + 
models. In Ref. ||l| it was shown that the cross section for sneutrino production pp — > Ue^e is 
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comparable to clgl production, and ei^e production is larger by a factor of 2-3 for a fixed 
value of mg^. Thus, the viability of the 6677 + event as ez, production (and the ability to 
distinguish from cr) depends in part on the phenomenology associated with sneutrinos. 

The dominant branching fraction of sneutrinos depends on the size of the Z component 
of the neutralinos and the gaugino mixings of the chargino, in addition to the mass hierarchy. 
There are 4x3 kinematic possibilities, where rrii, is lighter or heavier than "i^^^^ and rrifj^^. 
In the limit of pure neutralino states A^i = and = 7; the sneutrino has no coupling to the 
lightest two neutralinos since it does not couple to either pure state. Thus, in the case where 
mi, < m^^ ^ , the dominant decay of v will be to the kinematically accessible neutralino with the 
largest Z component. The relative branching fraction into Ni or N2 is therefore determined by 
the size of their Z component impurity. The branching ratios are shown in Table ^. 

For ei models, mj> > m^^, so that decays into the lightest chargino are always possible. 
The branching ratio for v — > Cil is always larger than 53%, while the branching ratio for 
the D — > C2I channel (if open) can reach 26%. The next largest channel is — > Niv, with a 
branching ratio up to 36%. The decay v N2I' is always open, but with a branching ratio 
below 6% due to the small Z component in 

For e/j models mj; is unconstrained, so the decay v ^ Niv \s the only mode that is always 
open. If decays into A^2 are also allowed, then the dominant decay of v can be into either Niu 
or N2i'. In special cases, we found it is possible for the Z impurity to be larger in than 
A'^i, thus the dominant decay could be P — > N2V. This is possible when < m^^ < rn^^, i.e. 
decays into charginos must be kinematically forbidden (an impossible scenario in cl models). 
When a channel into a chargino is sufficiently open, it dominates over decays into the lightest 
two neutralinos by a factor of more than 10. However, if the sneutrino is heavy m/^ > m^y^^, 
decays into the heavier neutralinos can be moderately large (branching ratio 10-30%), with 
decays u — > N^v dominating over u N-^v due to the larger Z component in A'4. 

4.6 Selectron branching ratios 

We have already discussed cl branching ratios for ei models, and branching ratios for 
cr models in Sec. since they are a fundamental part of the model building. The other 
slepton [en in cl models, and sl in models), will have branching ratios similar to ei (in ei 
models), or cr (in eR_ models) if its mass is roughly included in the 6677 + allowed range. 
In general, cr, cl will decay into the kinematically allowed final states with neutralinos, with 
the largest branching ratio for the channel N2e, if open, ei can also decay into Ci^2i^e if open, 
with a maximum branching ratio of 27% and 59% respectively. 
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Final State Kinematic Condition Range in cl models Range in e/j models 



NiUe - ^36 ^ 100 

N2i'e rriy^ > rrifj^ 5.5 97 

Nz^e rriy^ > "^7V3 ^29 ^22 

Cie mc^ > m^^ 53 ^ 94 100 

(726 m,>^ > -^26 ^48 



Table 9: Ranges of selected branching ratios (in %), as in Table |. 
4.7 Predictions for LEP 

The imminent upgrade of LEP to ^/s = 161 GeV (LEP161) and the forthcoming upgrade to 
^/s = 190 GeV (LEP190) provide a potential testing ground for the models constructed. With 
expected integrated luminosities of 25 pb~^ and 500 pb~^ (per detector), the one event level is 
at 40 fb and 2 fb for -y/s = 161, 190 GeV respectively. The first priority is to identify which 
processes have non-negligible production cross sections, then determine the possible signatures 
that depend on the branching ratios of the produced sparticles. It is important to emphasize 
that the following predictions assume the minimum cut ^ = 5 fb is placed on the a x B'^ 
for the 6677 + event. For instance, in some cases we are able to predict a non-negligible 
minimum number of events with a particular signature must be produced, although we do 
not necessarily give detector efficiencies. In principle, if one could demonstrate that failure to 
detect such events implies they do not occur at all, then only two possibilities remain: (1) A 
supersymmetric explanation of the 6677 -|- event in our framework must rely on an upward 
fluctuation from a x even lower than 5 fb, or (2) a supersymmetric explanation in our 
framework is not possible. 

Based on Observation 4 in Sec. |^, selectron production is always kinematically forbidden at 
LEP161 and LEP190 for the selectron that satisfies the kinematics. The other slepton (ept in ei 
models, or cl in models) can potentially be kinematically accessible at LEP161 or LEP190 
by simply requiring its mass be less than the threshold. This is obviously not a requirement (nor 
a constraint) of the selectron interpretation of the 6677 -|- event, and so we ignore selectron 
production at LEP. However, in e^, models it was shown in Eq. (|6|) that mj^ must be less than 
mg^, and so sneutrino production could be a visible signal at LEP190 (since rrij}^ > 81 GeV for 
all cl models), as will be discussed below. 
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4.7.1 LEP161 



In Fig. 11 we present all of the chargino/neutralino production processes that have cross 
sections above about 10 fb. The cross sections were computed with initial state radiation 
effects included. In ei or ei + sr models, none of the processes need to have large cross 
sections, although if it were possible to establish an upper bound on a^NiN^) < 600 fb, then 
an upper bound onaxB^ for the ee^^ + ^T event can be established at 7.5 fb, and in cr models 
o" X < 5 fb (i.e. all of our cr models would be excluded). Given the cut ^ = 5 fb, then in cr 
models one expects a minimum of 22 NiN^ pairs to be produced, but no other process (nor any 
processes in or ei + sr models) can have non-negligible minimum rates at LEP161. There 
are only four processes that could have large rates, which have the following maximum 

e+e- NiNs (55, 56, 49) 

N2N2 (19, 22, 12) 

N2N3 (11, 16, 7) ^ ' 

CtC{ (48, 132, 42) pairs produced 

for (e^, cr, cl + cr) models. Notice that the maximum pair production rates are always largest 
for cr models, then models, then + cr models. The rate for NiN^ is roughly the same 
in all models since the cross section is dominated by Z exchange. For the other processes, 
differing interference contributions between the Z exchange and light slepton exchange cause 
the differences in the production cross sections (see Table In addition, stop pairs could be 
produced at LEP161 (see Table |l] 



The character of the signal from N1N3 production is completely dependent on the decay 
of N3 which was described in Sec. (see also Table The dominant decay possibilities are 
N3 iVi"Z", N3 NiA{h) (if mAirrih) < 60 GeV), N3 IrI in ei models (if m^-^ < m^J, 
and N-i vv in eR models (if mi, < ^7^^^). The general signature is therefore Extra 
hb + ]fi occurs if the mass difference rrif^^—mfj^ is larger than rrih or rriA- Some other signatures 
are possible in special cases: In ei models one could have excess I'^l^ + ]p (if m-^^ < rriff^), or 
'yl~^l~ + (if mj^ < m^^). In cr models the decay — > i/i^ becomes the dominant decay if 
the sneutrino (and necessarily II) are light. Thus the dominant signature could be invisible, 
or 7 + ^, or + ^, if the mass hierarchy is mc < {ruff^, rriQ^, or < mt> < m^_^, or 

< my respectively. However, in these cases the cross section for IlD at the Tevatron would 
be quite large (see below). 

The dominant signal of N2N2 production is 77 + in all models. Note that the process 
a{N2N2) is always accompanied by a{N2N3) at a comparable rate (when kinematically allowed), 
which has the same signatures as N1N3 production (as above) plus one photon. 

CiCi production can be present with a large rate, the decay signature of Ci being the 
usual ^W"+^ in all models (see Table |^. The exception is if (and possibly rn^^) is 
lighter than m^_^, which can happen only in cr models. In this special cr model scenario, if 
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Figure 11: The range of the non-negligible cross sections at LEP161, for all models (shaded 
bar on left) and only br models (thick solid line on right). Each bar represents a particular 
production cross section, where the maximum and minimum height of the bar (or thick solid 
line) is the maximum and minimum cross section respectively. The shading on the left bars 
indicates the range of cross section for all models passing the cut A = 5, 7.5, 10 fb. Bars that 
touch the X-axis correspond to cross sections that can be smaller than 1 fb. 
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nil, < '"^c^) then the decay signature is hkely invisible. However, if the decay P — > N2I' is large, 
then the signature is 77 + ^. If m^^ < m^^, then additional possible signatures are + ^ 
(if m^-^ < rn^^), or + ^ (if m^-^ > rrij^^). Notice that the latter could be an additional 

source of 6677 events (see Appendix^). These remarks assume the stop is heavier than Ci. 

As an aside, we find that a maximum of (14, 13, 12) NiNi pairs can be produced, which 
can be observed as a 7 + signal once visible initial state radiation is attached. Although the 
SM background is severe, there are other possibly important contributions from e.g. jNiN^i^-^ 
vVNi). 



In Fig. 12 we present the ranges of the inclusive production of particular signals at LEP161 
for ei and e/j models. These signatures were generated by searching all possible decay paths. 
No efficiencies resulting from detector geometry or lepton/photon energy cuts are included. If 
the signals are the result of decays with moderate mass splittings, then presumably some of the 
events could be detected after applying reasonable cuts. A lepton I can be either e, or r, with 
either charge ±1. In particular, when referring to a "2/" signal, we sum over all family and 
charge possibilities (including, e.g., like-sign dileptons). X can be any combination of leptons, 
photons, jets, or nothing. In addition, all the signals implicitly include missing energy in their 
signature. We only include chargino/neutralino production processes in the inclusive sum, since 
eL in ei models and en in cr models is too heavy to be produced. If the other slepton {en in 
models, in models) is light, then the maximum cross section for particular signatures 
can be higher. 

Jet production is also an important signal. If N1N3 production is kinematically allowed 
and if only 3-body decays of N3 occur, then the rate into the jj + X + ^ signal is between 
roughly 400-1800 fb for both ex, and models. If chargino production is open then the rate 
can be larger. But, if 2-body decays into sleptons are open for A^3, then the rate can be near 
zero. 

Notice that in cr models only the 21 + X{-\-^) must be produced, the rate being between 
~ 2 to 20 events in 25 pb~^ of integrated luminosity. The reason that the 21 rate always 
has a non-negligible minimum is due to a combination of effects: cr models have a minimum 
a{NiN3) > 850 fb, and decays of Nil^l~ are always non-zero, even if 2-body decays 

operate. If only 3-body decays occur, then N3 Ni^^Z"{^ occur, with a rate of nearly 

10% (summed over families). Alternatively, if N3 vu is open, then u — > IC{-^ Iv) is the decay 
pattern. If rriu < m^, then it turns out that m^^ < m^^, and so decays N3 1^1,2) 
are non-zero, giving an appreciable 21 signal. All of the other inclusive signals could have rates 
smaller than the one event level. If one of these signatures were found (and deduced to be 
above background), then looking in the other channels might serve to confirm the signal. 

One promising signal is 77 + without any other event activity, which primarily originates 
from N2N2 production in the selectron interpretation. (This is part of the inclusive signal 
77 + X + ^ described above.) In a scenario with a gravitino LSP, we found that the standard 
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Figure 12: Range of inclusive cross sections for selected signatures at LEP161 without detection 
efficiencies; all signatures necessarily have missing energy in addition to that above. The shaded 
bar on the left corresponds to cl models and the thick solid outline on the right corresponds 
to en models. Here, X = leptons, photons, jets, or nothing, and Z = e, /x, or r summed over 
both charges and all three families. See the text for details. 
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Missing Invariant Mass (GeV) 

Figure 13: Comparison of the missing invariant mass distribution in the 77 + signal at 
LEP161 from two different selectron interpretation models, (a) a sample A'^i =LSP model with 
37,65 GeV (dashed line), and (b) a model with a gravitino LSP and =65 GeV 



m 



(solid line). 



model background for 77 + ^ is distinguishable from the gravitino signal e^e^ NiNi 
jjCrG using the missing invariant mass distribution [Q]. Here, we point out that a selectron 
interpretation with a neutralino LSP can be distinguished from one with a gravitino LSP using 
the missing invariant mass distribution, assuming the SM background is small (see Ref. [Q] for 



a discussion of the background). In Fig. 13 we show the missing invariant mass distribution 



Md 



{Pe+ + Pe- ~ P'yi "^72) LEP161 for two different models: (a) The cl sample model 



37, 65 GeV, and as usual A'^i 



i i^b, — 7- (b) A model with 
7. The difference in the missing 



in Appendix ^ with m 

a (very light) gravitino LSP with =65 GeV, and A^; 
invariant mass distribution illustrates how the scenarios might be distinguished using the 77+-^ 
signal. It should be noted the general character of the missing invariant mass distribution for 
the gravitino LSP model in Fig. 13 is not particularly sensitive to nij^^, but simply that is 
very small compared to the neutralino or selectron masses. 



4.7.2 LEP190 



In Fig. 14, we present all the chargino/neutralino production processes with cross sections 
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possibly larger than about 1 fb for LEP with y/s = 190 GeV. As above, the cross sections were 
computed with initial state radiation effects included. Now, N1N3 production must be large in 
all models, and many other processes can easily give large rates. The processes with large rates 
include all of the ones at LEP161, and also N1N2, NiNa, N2N4, CfC^. The maximum rates 
are as follows: 

e+e- ^ N1N2 (20, 6, 24) 

iViiVg (785, 780, 780) 
iViiV4 (82, 79, 78) 
N2N2 (505, 560, 346) 

iVaiVg (335, 416, 230) ^ ' 

N2NA (73, 64, 34) 

C^+Cf (965, 2120, 1195) 

CtC^ (409, 695, 350) pairs produced, 

for (cl, Br, cl + gr) models. For N1N3 production, the minimum number of pairs produced 
is (400, 475, 320) for (e^, eji, ei + eji) models given the minimum threshold ^ = 5 fb. For eji 
models only, a minimum of 5 NiN^ pairs, 25 N2N2 pairs, 40 -/V2-/V3 pairs, and 250 CiCi pairs 
must be produced given the minimum threshold ^ = 5 fb. As for NiNi pair production we 
found a maximum of (177, 164, 152) pairs can be produced. 

The detection signatures for the chargino/neutralino pairs common to LEP161 are the 
same as above. Here we discuss the processes that are different. First, the process N1N2 
gives a 7 + ^ signature. The signatures for NiN^ and N2NA are entirely dependent on the 
Na branching ratio; Na can decay in a variety of ways outlined in Sec. |4.3| . Perhaps the most 
striking signature is when A^4 Nih{A), giving a bb + ^ signature for NiNa production and 
667 + ^ signature for N2NA production. The signature of the process CfC2 also depends 
crucially on the branching ratio of C2, but one lepton with perhaps one photon plus missing 
energy is typical (assuming the stop is heavier than Ci). Thus, a reasonable expectation for 
C^C^ production is + fJ. It is also possible that only 3-body decays of C2 are open, 

in which case no photon would appear in the final state. The final states from C2 decay are 
summarized in Table |^. 

In addition, sneutrino pair production (if open) is another process that is relevant for 
cl models. To have D^De production kinematically accessible with nic,^ < 95 GeV, then the 
sum rule in Eq. (^) implies tan/3 > 1.2 is required for mg^ > 100 GeV (as needed by the 
kinematics of the 6677 + event), and for mg^ = 107(118) GeV, then tan/3 > 1.5(2.8). 
Hence sneutrino production in ei models never occurs at LEP190 if mg^ > 118 GeV. The 
signature of Uei^e production depends on the sneutrino branching ratio, but it was already 
established in Sec. that Ue eC'i(— *■ iVi'W) is the dominant decay pattern. Thus the 
signature is ee'W 'W'+^, which is indeed quite prominent. 



In Fig. 15 we present the ranges of the inclusive production of particular signals at LEP190 



for cl and e/j models. As in Fig. no detection efficiencies are included. Notice that while only 
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NiN^ production had a non-negligible minimum rate (see Fig. |TJ), both the signals 21 + X + ^ 
and 77 + X + ^ (which rarely comes from N1N3 production) are always larger than one event. 
Further, inclusive production of /7 + X + ^ and llj + X + ^ are always larger than the 10 event 
level for cr models only. All of the other inclusive signals could have rates smaller than the 
one event level. As in LEP161, if one of these signatures were found (and deduced to be above 
background), then looking in the other channels might serve to confirm the signal. Another 
important search strategy would be inclusive signatures with jets (+ photon(s)) that can have 
significantly larger rates than the lepton(s) (+ photon(s)) signatures. 



4.8 Predictions for Tevatron 



The assumption underlying the selectron interpretation is that the Tevatron has already 
observed a candidate selectron pair production event. Because many more states of the under- 
lying supersymmetric model are accessible at a hadron collider, we here focus on the associated 
signals that should be observed in the present data set (100 pb^^ per detector) or in the next 
scheduled upgrade (1-2 fb^^ per detector). As in Sec. 4.7, we identify the processes that have 
non- negligible production cross sections, then determine the possible signatures that depend on 
the branching ratios. Again, it is important to emphasize that the following predictions assume 
the minimum cut ^ = 5 fb is placed on the a x for the 6677 + event. 

In Fig. ^ we present all of the chargino/neutralino production processes that can have cross 
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Figure 15: Range of inclusive cross sections for selected signatures without detection efficiencies, 



as in Fig. 12, but for LEP190. 



sections above about 50 fb. We use leading order CTEQ3L structure functions evaluated 
at = s. At the Tevatron the cross sections do have a contribution from t-channel squark 
exchange (see Table |2|), but the dependence on the squark mass is usually weak for the squark 
masses in our models. If only eji is light, then these are the only necessary associated processes 
to the 6677 event. If is light, however, then there must also be eL^e and Pg^'e production. 
In Ref. Q we found cj(t'et'e) ~ cr(eLeL) and a{eLi^e) ~ (2 — > 3)(T(eLei,) for the same mg^, i.e. 
the cross sections are typically tens of fb. It is also possible that both and cr can be light; 
in particular, the other slepton (en in cl models, bl in br models) can be lighter than the one 
giving the 6677 + event, which can dramatically affect the signatures. The pair production 
processes that have the largest cross sections and also have a non-negligible minimum cross 
section are given in Table |l^, where the full range from the minimum to the maximum number 
of pairs produced for an integrated luminosity of 100 pb^^ are shown. 



The signatures for NiN^ and C^C^ are the same as for LEP (described in Sec. |4.7.1 ) 



and the decays of C2 were also discussed in Sec. 4.7.2| . For completeness we list the possible 



signatures of all of these processes here: will mainly give or bb + if "m^A < 

60 GeV. If 2-body slepton decays are allowed, then in bl models one can have + fJ^, or in 
Br models one of invisible, 7 + ^t, or l^l^ + ^t- C^Ci production gives typically + ^7-, 
or if 2-body decays into /l, occur (in br models only), then depending on the mass hierarchy 
one can have 77 + -^^; or + ^t, or 1~^1~^j + ^t- C2C2 production gives similar signatures 
as CfCi production, given — > m^j^ and allowing for the possibility of 2-body decays in 
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Process Range in models Range in models Range in ei + gr models 





31 - 


129 


43 - 


^ 145 


29 - 


^ 128 




40 - 


285 


56 - 


^ 264 


29 - 


258 




8 - 


85 


28 - 


79 


15 - 


77 




75 - 


638 


132 - 


540 


54 - 


552 


C^N2 


2 - 


75 


3 - 


80 


1 - 


75 


CtN3 


32 - 


98 


36 - 


103 


28 - 


96 


CtN2 


2 - 


76 


15 - 


69 


5 - 


74 




3 - 


51 


17 - 


54 


8 - 


55 



Table 10: The range of the number of chargino/neutralino pairs produced at the Tevatron 
assuming an integrated luminosity of 100 pb^^. The processes displayed here include those 
that have both a large production rate and a non-negligible minimum production rate. 

the context of both cl and en models as above. 

The processes C^Nj are unique to the Tevatron, with CfNi, C^N,, C^^s^ CtN2 and 
giving the largest rates. As described above, the chargino typically gives jj + IpT and 
+ fT, although possible 2-body decays into sleptons can give 7 + or I^ + ^^t, or I^'j + ^t- 
Thus the signature of C^Ni production is one of the above signatures for a single chargino. The 
signatures of C^N, and C'^iV2 are as above plus one photon. Finally, the decays Cj^iVa and 
are one of the above signatures coupled with or A'4 decay. Here again we can utilize 
Sees. [1.7| , to obtain the possible decay signatures. For A'^a, the decay signature is 
bb + (if f^A < 60 GeV), and if 2-body decays to sleptons are open then for ei models the 
signature could be or 'yl~^l~ +$Ti while for models the signature could be invisible, 

or or +fT- Thus, if only 3-body decays were open for charginos and neutralinos the 

signature of CfN^^ and C^N^ would be 'W"'Z"+^t, which gives the well-studied trilepton 
signal [^]. If 2-body decays of the charginos or heavier neutralinos are present, then one or 
more photons could be present in the final state, with possibly fewer leptons. 



In Fig. 17 we present the cross section for many promising signatures at the Tevatron. As 



in Figs. 12 and 15, no detection efficiencies have been included. We include chargino/neutralino 
processes in the sum, as well as enen production in models, and clcl, t'e^'e) ^Li'e production 
in models. We see that all six inclusive signatures involving leptons or photons are expected 
to have minimum rates of roughly 2 to 30 events, regardless of the type of model {ei or e/j). The 
77+ A, Z7+A and /77+A signatures can be much larger in cr models, but this only happens in 
the particular kinematic scenario with < mj^^rriii < m^^(< rrii^). In this case, charginos 
always decay through the 2-body channels C — > Ili' and C ul, with liji^' — > N2{-^ A17). 
Thus, processes with intrinsically large cross sections such as CiNi production can lead to a 
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large Ij + fix signal, and similarly for other processes involving charginos. 

The Ij + ^T (and jjj + ^T) signals are important and can arise from: C1N2 and C2N2 
production in models with C — > IvNi; CiNi production in models with Ci — > IN2{-^ 
Ni'y)) or Ci — > /^(^ vN2{-^ Ni^)); and /^P production with 1^ — > IN2{-^ -^17)- The chargino 
decays assume m^^ > m^^. For just CiN2 production there are roughly 10-130 pairs produced 
in the present CDF and DO samples (each) with the probable signatures 7+'W'+^ (before 
cuts); 'W decays to jj or l^u as usual. For 'W'^ jj, these events have no parton-level SM 
background. 

Many of these signatures should be detectable, since the mass differences between super- 
partners is often constrained to be small but non-zero, as in Fig. |8|. For example, in decays such 
as iVs iVi"Z" and 6*1,2 ^ iVi'W", the invariant mass of the virtual "Z" or 'W can be large. 
In particular, the invariant mass of the "Z" from decay is between to 40-60 GeV, thus 
an excess in pairs of leptons (or jets) that reconstruct to an invariant mass mi+i- < 60 GeV 
accompanied by a large missing energy is a distinctive signature of NiN^ production in our 
models. 

In addition to classifying the most promising signatures, we have also performed a number 
of event level simulations for a limited subset of our ei and en models, with the other slepton 
heavy. The purpose is to get a feeling for the efficiency of detecting multi-lepton and/or 
photon signatures. First, we address the issue of efficiencies for the 6677 -|- event, since this 
is important for interpreting the threshold A in the 6677 rate. An efficiency represents the 
probability that a certain class of events passes a particular set of cuts defined before the data 
is analyzed. We chose a set of cuts such that: (1) the event would be triggered on and analyzed, 
and (2) the event would not suffer from obvious detector backgrounds like jets faking leptons 
or photons. To show the dependence of our efficiencies on the particular set of cuts, we choose 
a loose set with < 2, < 1, {p'^'^\$t) > EJ^'^ = 10 GeV and a tight set identical to the 
loose set except i?™™ = 20 GeV. The efficiencies we found range from 0.02-0.23 for the loose 
cuts, and from 0.01-0.12 for the tight cuts, but efficiencies outside these ranges (from models 
not covered in the subset) are possible. If E'™™ is increased to 25 GeV, the mean efficiency is 
0.04. The loose cuts are sufficient for CDF to have triggered on the 6677 event. 

We have also studied Ij, II, and llj signatures using a similar set of cuts (^r > 20 GeV and 
pj'^^ > Ef""), where for the purposes of detection I is summed over e and fi only. Typically, 
when = 10 GeV one expects between 1-5 (2-12) I'y + ^x events in 100 pb^^ for Br models 
{cl models) from chargino/neutralino production alone. An additional 1-2 events are expected 
from blDe production in ei models. This result is essentially unchanged for the simulation subset 
of models if i?™™ = 20 GeV. This is expected at least for the photons since the kinematics 
enforce hard photons in the final state from slepton decay. The SM background from W'j 
production yields 105 and 37 events for each set of cuts respectively. With tighter cuts, it is 
possible to achieve a signal to background ratio near one for some models. The expected II 



41 



CO h-H 

O TO 



05 



o >> 

EC M 

Qrq 
o ■ 



cross section [fb] 



cross section [fb] 



o 

5' 





cm 



< 

1-i 

o 

1^ 



o 




00 

o ^ 



H- 
P' 



C1N3 



C1N4 



C2N1 



C2N2 



C2N3 



C2N4 



N1N2 
N1N3 



N1N4 



N2N3 



C1C1 



C2C2 



1-i 

o 
o 

fB 
CD 
CD 
CD 
CD 



10^ 



^ 10 

g 
'■♦-« 
o 

0) 
CO 

w 

CO 

o 



O 



10' 



10' 




X 


X 


X 


X 


X 


X 


+ 


+ 


+ 


+ 


+ 


+ 


cvi 










CO 















Figure 17: Range of inclusive cross sections for selected signatures without detection efficiencies, 



as in Fig. 12, but for Tevatron. 



signature, resulting mainly from CiCj production, is between 0-6 events for Ej^™ = 10 GeV, 
and 0-2 for E"™™ = 20 GeV. Similarly, the expected ll'j signature is between 0~5 events for 
E^^^ = 10 GeV, and 0-2 events for = 20 GeV. Other signatures, such as 77, ^77 and 

31, produce at most 1 or 2 events for £'™™ = 10 or 20 GeV. Therefore, it would appear that 
the Ij channel is the most promising for confirming the supersymmetric interpretation of the 
6677 event (assuming m^^ > m^), though other signals with limited backgrounds are clearly 
possible. 



4.9 Alternative interpretation 



Throughout this section we have described the constraints and predictions in the selectron 
interpretation. However, in Sec. |2|we described an alternative interpretation involving chargino 
production that could explain the 6677 + event. Those readers interested in the model 
building associated with the chargino interpretation are referred to Appendix which provides 
many details and an example model. 
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5 Comments on models with a light stop 



We have seen that the effect of requiring a large a x for the 6677 + event is to 
strongly constrain the chargino, neutralino and slepton sections. Up to now, we have assumed 
the squarks are sufficiently heavy so as not to directly interfere with the necessary decay chain. 
However, it is possible that a light stop ti can exist simultaneously with the needed hierarchy 
in the other sectors. In particular, neutralino decays — > tit are absent in our models (with 
the cut ^ = 5 fb), since all neutralinos are lighter than the top quark. Therefore, the decay 
chain in the selectron interpretation need not be disrupted, if the radiative neutralino decay 
can be large with a light stop. 

The chargino interpretation described in Appendix ^ is a different matter, since charginos 
would always decay to the light stop C — > ti6 if kinematically accessible. This is true regardless 
of the mixing angle 9f: that determines the W^-i coupling, since the Yukawa coupling H^- 
t is large. Thus, it would seem that a chargino interpretation of the 6677 + event from 
pp C^Cj is not possible unless m^^ > m^. This is basically the scenario described in 
Appendix 

To construct models with a large 6677 + event rate and a light stop, one must consider 
the effects of a small m^^ on the radiative neutralino decay width and on the mass hierarchy. 



As we have remarked in Sec. 4.2, the dynamical mechanism for a large radiative neutralino 
branching ratio appears not to be strongly dependent on m^^ [^]. For instance, models can 
be constructed with m^^ = 50 GeV, m^^ > 250 GeV, and a large radiative neutralino branch- 
ing ratio arising from the dynamical mechanism. However, some suppression to the radiative 
neutralino branching ratio from light stops is present, so the 6677 + IpT rate is maximized 
in the limit of all squark masses large. For example, the largest 6677 + Ifix rate in cl and 
models with a light stop is 13.8 fb and 6.1 fb respectively. Since < m^-^ must be 

obeyed so that A^i =LSP, the upper limit on m^^ can be more restrictive than found above if 
nii^ < 74 (50) GeV in ei (en) models by Observation 1 in Sec. ^. This induces a rough upper 
limit on which also has implications for the chargino masses. 

There is an additional degree of freedom in the value of 9^, which determines the SU(2) 
couplings of ti^2 with the gaugino components of charginos and neutralinos. Maintaining a 
hierarchy between m^^ <C (— rriq) would seem difficult without giving a large 6p but 
this can be avoided if ti ~ t^ (or 0^ ~ tt/2 by our definition). However, requiring m^^ ~ 50 GeV, 
0^ ~ tt/2 implies m^^ must be large (perhaps of order 1 TeV or more) if the light Higgs h is to 
have a mass that is not excluded by LEP. In general this implies that rrih will lie within the 
region accessible to LEP, though further analysis is needed to be precise; sin^(/3 — a) can be 
below one, and nih can be near its present lower limit from LEPl. Note also % slightly affects 
the radiative neutralino decay [11|. 
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Another constraint on models with a light stop comes about if + rrib < mt. Then top 
quarks must decay into stops with a branching fraction of about 1/2 if ~ 50 GeV. It was 
observed in Ref. [0] that a branching ratio of t ^ ti6 of 50% is not excluded by Tevatron data, if 
gluinos and squarks with masses of roughly 0(250) GeV exist, giving additional top production 
to supplement the SM contribution while half the top quarks decay into the lightest stop. For 
our purposes we note that if the masses of non-stop squarks are greater than roughly 250 GeV, 
then they are not crucial in maintaining a large radiative neutralino branching fraction. 

The simultaneous existence of a light right stop, a heavy left stop, other squarks (except 6l) 
and the gluino with masses 0(250) GeV, and a large 6677 + rate is therefore an interesting 
possibility. We explicitly constructed nearly 200 models, mostly of the ei class due to their 
larger cross section. We did not find significant differences in the models' distribution in M1-M2 
plane, nor in the ^-tan (3 plane. However, regions in these planes that were populated by heavy 
stop models with a x near the ^ = 5 fb cut are no longer are allowed. For instance, no light 
stop models approached the gaugino mass unification (M2 = 2Mi) line, and was restricted 
to be less than 62 GeV. Hence, there are a number of phenomenological consequences of 
assuming a light stop (m^^ = 50 GeV). First, as noted above, the branching ratio of iib 
is virtually 100% (when kinematically accessible), followed by the one loop decay ti cNi 



if < rrifj. Thus all the signatures as noted in Sees. 4.7 and 4.8 arising from charginos 
become bc + ^T- For example, while the dilepton signal from N1N3 is unchanged, the dilepton 
sig nal from CfCj becomes hbcc + ]pT- Also, the restriction rrij^^ < m^^(= 50 GeV) results in 
somewhat tighter restrictions on the upper bounds of the other chargino and neutralino masses. 
In particular, ^ 90 GeV and ^ 100 GeV in our light stop models, and the sum 

{rrijy^ + rrij^^) ^ 150 GeV. One consequence is that N1N3 production is always kinematically 
allowed at LEP161, with a cross section in the range 1.1 < cj(A^iA^3) < 2.1 pb. 

Stop production at LEP161 may be directly visible with the expected integrated luminosity 
if 2m^^ is below threshold [28|. In Table ^ we present the cross section for stop production 



at LEP161 and LEP190 for a selection of light stop masses. At LEP161 one would expect 
roughly 20 (5) stop pairs produced per detector, for m^^ = 50 (70) GeV. At LEP190 one would 
expect roughly 380 (95) stop pairs produced per detector, for m^^ = 50 (80) GeV. All of the 
cross sections were calculated with approximate final state QCD corrections and QED initial 
state radiation effects included, and assuming ti = tR. Also, titj bound state effects can be 
important close to the threshold. 



It has been noted [16| that when there is a light stop (so that Cj — > ti6 and t tiNi), there 
is a large set of events predicted at the Tevatron by supersymmetry that has no parton-level 
SM background. Even after all branching ratios and detection efficiencies are included, tens of 
events remain in the present 100 pb~^ at Tevatron. These events arise from three sources, (1) 
Cf{^ iib)N2, see Table P; (fi) Wb)t{^ tiiVz); (ifi) qN2)q{^ qg{^ Wb)ii))- 

In all cases, ^i7; ~^ cA'^i, and typically W jj. After branching ratios and cuts 

there should be approximately 35-100 events with the signature 67 + ^+ jets. 'Jets' means 
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m^^^j^ Cross section (in pb) 
(GeV) LEP161 LEP190 



50 


0.85 


0.76 


60 


0.50 


0.56 


70 


0.20 


0.37 


80 




0.19 



Table 11: Cross sections for light stop ti (= tji) production at LEP161 and LEP190 with 
approximate final state QCD corrections and QED initial state radiation effects included. Close 
to the threshold the cross section values may receive large corrections due to titj bound state 
effects. 

1-5 parton level jets, including 1-2 charm jets (an average of 1.5/event). This prediction could 
lead to a sample that allowed a robust (rather than one event level) detection of superpartners 
in the present CDF and DO data. When W ^ Iv for these events, additional good signatures 
arise and one expects an excess of 'W6c events that would appear in the top sample, and 
+ ]p+ jets events. 

The simultaneous existence of a light stop and a light chargino (as necessarily arises in 



6677 + models) can give rise to a shift in i?5 [29]. We have analyzed models with nii^ = 
50 GeV, ti = tfi and find that the maximum shift in Rf, is 5Rf^^^ ^ 0.003 from chargino- 
stop loops only. Charged Higgs-top loops can also be significant, with a shift SRb < —0.0005 
depending on rriA- In all cases tan/3 must be near 1 for a maximal shift in Ri,. For example, 
tan /3 = 1.1, 1.5, 2.0 can all give a large 6677 + f^T rate, while the shift in R^ is at best 0.0028, 
0.0021, 0.0018 for chargino-stop loops only. Further, Ri, is sensitive more to the parameter tan f3 
than m^^ , as is clear since the chargino mass is inversely related to tan /?; in the above three 
cases m^_^ is roughly 83, 80, 70 GeV. We note that these calculations have been done assuming 
m^^ = 50 GeV, 6^ = tt/2, which is nearly optimal since the maximum shift in Rf, decreases as 
either the stop mass is increased, or 9^ is taken far from 7r/2. 

As has been emphasized, getting a significant shift in R^ requires a chargino that has a 
large Higgsino component, and the related result that fi is small and negative. It is interesting 
that the value of /U and the chargino properties coming from the analysis of the 6677 + 
event have the properties needed to give such an effect. Finally, we note that a shift in Rf, 
necessarily implies a shift in as extracted from the LEP Z lineshape, through the relation 
6as{Mz) ~ —4:6Rh This limits the maximum shift in Rf, to about 0.0025, consistent with 
the above numbers and giving Rf, < 0.2182. It is worth emphasizing that a significant shift 
in Rfy (and as) is only possible simultaneously with a supersymmetric interpretation of the 
6677 + event if Ni is the LSP [|]. 
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6 Concluding remarks 



We have seen that supersymmetry with Ni =LSP is a viable explanation of the CDF 
6677 + event. The primary constraints are the kinematics of the 6677 + event, the 
radiative neutralino branching ratio ^^il, the selectron decay e — > eN2 and LEPl— 

LEP130 data. Given a minimum threshold on the cross section times branching ratio of pp 
e^e~ e~ N2N2 — > e'^e~^^NiNi at the Tevatron, a selectron interpretation requires Mi, 

M2, /Li, tan/3, me in tight ranges (see Table ^ and Figs. |^, |8|). The corresponding chargino and 
neutralino masses and the cross sections at LEP and Tevatron are similarity constrained. This 
is the origin of the predictions made for both LEP and Tevatron based solely on the 6677 + ]pT 
event, where many signals can be large, and some must be produced. These signals are deduced 
from the cross sections and branching ratios without efficiencies, although in many cases the 
mass differences between sparticles cannot be arbitrarily small, and so presumably the signals 
are detectable. For example, N1N3 production must occur at LEP190 with the mass difference 
40 < — < 60 GeV in all models, which implies a pair of leptons or jets from the decay 
N3 Niff would have an invariant mass up to roughly 60 GeV. The inclusive signals that 
must be produced at LEP190 with an integrated luminosity of 500 pb~^ are a{2l+X + ^) > 50 
events, and o"(77 + X + ^) > 3 events. At the Tevatron, the inclusive signals that should 
have been produced (with an integrated luminosity of 100 pb~^) are a{2l + X + Ipx) ^ 30 
events, £1(77 + X + $t) Z 2 events, a{l-f + X + $t) Z 15 events, cj(2/7 + X + $t) Z 4 
events, a{ljj + X + ^t) ^ 2 events, and a{3l + X + J^t) <L 2 events. All of these signals 
assume X = anything (leptons, photons, jets), and are valid for ei or e/j models. For cr models 
only, the inclusive signal that must be produced at LEP161 with an integrated luminosity of 
25 pb~^ is a{2l + X -\- ^) > 2 events. Also for eji models only, the inclusive signals that 
must be produced at LEP190 (in addition to the ones above) are 17(^7 + X + ^) > 5 events, 
and (7(2/7 + X + ^) > 5 events. We have examined many inclusive signals with leptons and 
photons, but of course inclusive signals with jets (+ photons) are also important and in some 
cases can be larger. 

The selectron interpretation can be made with the selectron e being either e^, e^j, or a sum 
over cl and cr contributions. The difference between ei and cr is in the SU(2)i couplings of 
cl, causing for example the cross section at the Tevatron a{pp e^eL) ^ 2.2a{pp — > e/je/j) 
(in the mass range of interest), and the presence of e^, couplings to charginos. Thus one way 
to distinguish ei (and ei + eji) models from cr models is with the associated charged current 
channel pp — > e^Pe that gives at least Ij + ^t, with possibly more leptons or photons depending 
on the decay of i)e- Studies of such signals are relevant for I = e, fi, t. A further source of 
l^^ + X + events comes from C'i^2(^ iVi/i/)iV2(— > Ni^), as well as CfCj with C ei^e or 
C — > z>ee if mg^ or mc^ is lighter than the chargino. Thus if no excess of associated events can 
be attributed to the absence of CL^e production, then it becomes less likely that the original 
selectron was e^, though it cannot be definitive until a clean result is published. In addition. 



47 



particular signals must be produced at LEP161/190 for eji models that are not necessarily 
present for ei models, and thus if these associated events were not found, then it becomes less 
likely that the original selectron was en, with the same caveat as above. At LEP it is necessary 
to study the relative rates of different channels to distinguish ei from cr, unless selectron pair 
production is actually observed there. In fact, if cl or eji production is observed (and the 
LSP can be established to be A^i), then we immediately know which charged slepton is not 
responsible for the 6677 + event, since as we have shown in this paper the slepton giving the 
6677 + is kinematically forbidden at LEP161 and LEP190. Thus there is no unique signal 
to discriminate from cr (from e/, + eji) models; only through the pattern of multiple signals 
can the nature of the selectron be determined. 

We have also seen that a chargino interpretation of the 6677 + fJx event is a distinct possi- 
bility. In either the selectron or chargino interpretation we expect at least the constraints from 
radiative neutralino decay to hold, and light sleptons are probably also a shared requirement 
for either interpretation (see Appendix^). One way to eventually distinguish the selectron in- 
terpretation from the chargino interpretation is to compare the rates of 6677, fifi'yy and 6/^77. 
Assuming a mass degeneracy among the sleptons of different families, the selectron interpreta- 
tion predicts roughly an equal number of 6677 and /i^77 events, with a significantly depleted 
e/i77 signal originating only from f~^f^ production followed by t^t~ — > e^/i^ -|- X. Alter- 
natively, in the chargino interpretation one would expect roughly double the number of 6/^77 
events as compared with either 6677 or ^/i77 events. Thus comparing the e;U77 rate with either 
6677 or A'/^TT would provide a useful means to discriminate between the two interpretations. 
Notice also that events of the type /+/'~77 -|- can be produced only from CfC^ and f+f" 
production. 

It is important to remark that the ee-^^+^x event phenomenology could be connected with 
other phenomena. If the LSP= A'^i is stable, then it could provide a cosmologically significant 
relic density even if it is mostly a Higgsino (as required by the 6677 -|- event). For a 
given value of Jl/i^ the mass of Ni is correlated with tan /?, and so gives a subset of the models 
constructed here. The predictions for associated phenomenology are tighter; and generally the 
signals can be larger. Also, we have described in detail the effect of assuming a light stop in 
addition to the 6677 -|- event, in particular its connection to (and as [@). A light 

stop has many other consequences m, that we will not go into detail about here. 

However, it is perhaps useful to remark on how can we learn if there is a light stop. The 
easiest way would be to observe it at LEP. The cross section ranges from about 0.2-0.8 pb 
over the range 50 < m^^ < 70 (80) GeV of most interest for LEP161 (LEP190). LEP190 with 
tens of pb~^ will be definitive. For such light stops and even for somewhat heavier ones up 
to nit — rrij^^ ( ^ 100 GeV in models considered here) searches in t{-^ tiNi)t{-^ Wb) or stop 
pair production can be definitive. Indirect evidence for a light stop before there is definitive 
collider data could come from a convincing Rfj excess, from slepton pair production at the 
Tevatron without associated leptons, and photons from chargino channels because Cf^ — > iib, 
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from anomalous behavior of top properties [^], and from the 7 + 6+ jets events [16| commented 
on in Sec. H- Note that ii could be near Ni in mass, and therefore give very soft fermions plus 
large missing energy. 

We have stated that certain signals must be produced at LEP and Tevatron, and some 
signals might be produced if kinematically accessible. For example, at LEP161 three neutralino 
and one chargino pair cross sections are large enough to give a signal if about 25 pb~^ is collected. 



The signatures are described in Sec. |4.7.1| and can sometimes be somewhat unusual. At LEP190 
many more processes can be open, which can all give signals with possibly unusual signatures 
(see Sec. 4.7.2| ). It is important to emphasize that the predictions assume the minimum cut 



^ = 5 fb is placed on the a x B'^ for the 6677 + event. In principle, if one could demonstrate 
that failure to detect the signals implies they do not occur at all, then only two possibilities 
remain: (1) A supersymmetric explanation of the 6677 + event in our framework must 
rely on an upward fluctuation from a x B'^ even lower than 5 fb, or (2) a supersymmetric 
explanation in our framework is not possible. We note that even if the cut ^ = 5 fb needs to 
be relaxed, there are still constraints from requiring a moderate branching ratio for — * A^i7 
as demonstrated in Fig. 

The 6677 + event has given us a profound example of how low energy supersymmetry 
could be discovered with one event. It is not obvious that such an event could be explained 
by supersymmetry, and we emphasize here the predictability of the theory once such an expla- 
nation is adopted. In particular, we have shown that assuming the 6677 + event is due to 
supersymmetry with a A'^i =LSP imposes strong constraints on the supersymmetric parameters, 
and predicts much associated phenomenology. Confirmation at LEP or Tevatron from the myr- 
iad of associated signals described in this paper is necessary to be definitive. It is remarkable 
how much can be learned from the Tevatron data, if the signal is confirmed. 
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Appendix A: Models in the chargino interpretation 



The chargino interpretation purports to explain the 6677 -|- event through chargino 
production and decay, a priori sharing only the requirement of radiative neutralino decay with 
the selectron interpretation. The possible sources of 6677 in the chargino interpretation are from 
pp — > C^CJ with Cij —>■ N2ei', followed by N2 -/V17. The decay Cij N2ev can proceed 
through either on-shell or off-shell W , ei and i^g- However, the 2-body decay C\ WN2 
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is not possible, since — m^_^ < 25 GeV when the radiative neutrahno branching ratio 
B{N2 -/V17) is required to be large. 



A.l Chargino production and 3-body decays 



If cl, z/g are heavy, the branching ratio for the decay Cj — > A^2efe is dominated by W- 
exchange, with a branching ratio the same as that for the SM decay W — > eug equal to 11%. 
Chargino production with heavy sleptons therefore implies for every event, roughly 20 

other events with jet activity (possibly accompanied by one charged lepton) or two different 
charged leptons (plus two photons). In addition, the channel Ci Niev is always open 
and it is generally favored by phase space, in particular in the case i = 1, since the mass 
difference m^^^ —mf^^ is never large. Further, it seems difficult to find a region of the parameter 
space allowed by LEP data, consistent with a large neutralino radiative decay branching ratio 
and the general kinematical 6677 requirements, where the non-radiative channels into A'^i are 
dynamically suppressed. This holds for both on- and off-shell VF-exchange, and as a result the 
branching ratio for the decay Ci — > N2ev hardly exceeds 6% for i = 1 and is even lower for i = 2. 
Hence, to get > 5 fb 6677 signal from CfCj production and decay Cij W'-*\^ eu)N2{^ 
7V17) one needs a cross section at least roughly 1.5 pb (even assuming B{N2 — > -/V17) = 100%), 
since >eZ^e) is well below 1%. This does not seem to be possible with 

an individual chargino pair production process, given all the other 6677 constraints. However, 
a small but non-zero signal can always arise from this source in models which are compatible 
with the selectron interpretation. We have found models with up to ~ 1 fb 6677 total signal 
from the sum of C^CJ production and 3-body Cf^ decay in our selectron interpretation models. 
(These contributions were not included in the selectron interpretation.) 



A. 2 Chargino production and 2-body decays 



We consider in the following chargino production followed by 2-body decays into sleptons, 
which allows an enhancement of the total possible branching ratio into the 6677 final state. The 
regions are somewhat different in the chargino interpretation with C ^ II than in the selectron 
interpretation; in particular we found the constraint raf;^^ — mfj_^ > 20 GeV is no longer 
required. (We have checked that a neutralino mass difference of order 10 GeV can be sufficient 
in the chargino interpretation.) This may in principle allow the kinematical mechanism for 
the enhancement of the radiative neutralino decay branching ratio to operate simultaneously 
with the dynamical mechanism to obtain a large 6677 rate. In Sec. 4.2 we already encountered 
particular models in the selectron interpretation where the kinematical mechanism plays an 
important role, and this may be true for the chargino interpretation to an even greater extent. 
However, a small mass difference w-^^ ~ ''^Ni seems only to be allowed when m^y^ is small, 
so that it can presumably receive a large boost after the or De decay and generate a hard 
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photon. The only way to construct a model with two very light neutralinos and a heavier 
chargino is to enter the "light gaugino-Higgsino window" (see Sec. 4.2), but even there it seems 



difficult to build a model which falls in the region suggested by the 6677 + event kinematics 
with the constraints from the branching ratios. Also, with a neutralino mass difference of order 
10 GeV or more the radiative neutralino branching ratio never approaches 100% from only the 
kinematical enhancement ||ll|. Hence, as in the selectron interpretation it would appear that 
the dynamical mechanism for a large radiative neutralino decay is required. This, along with 
the following argument for the need of a mostly gaugino A''2, explains why it seems possible 
to build models with large 6677 rates in the chargino interpretation only in regions of the 
gaugino-Higgsino parameter space similar to that in the selectron interpretation. 

The maximum B[C^ IN2)] for I = D^, is 1/3, and for I = ei is 1/6 due to the 

slepton mass degeneracy assumption and m/^^ < mg^ (assuming the decay into sneutrinos is 
not strongly suppressed). Also, the slepton decay channels with A^i in the final state are always 
open and enhanced by phase space. Thus to maximize the branching ratio into A''2, one has 
to minimize the N2 Higgsino components (which do not couple with sleptons) and maximize 
the Higgsino component of A^i. In this way, the branching ratio for bl is enhanced, 

analogous to the selectron interpretation. Typically, the branching ratio for the combined 
decay CfCj — > 6677 though 2-body decays into sleptons can reach at best ~ 4%, assuming 
B{N2 Ni'j) = 100%. In the C2 case, a further source of suppression can come from the 
channel Cie (if open), that always dominates over U(, N2i'e or z>e — > A'^ii^e- A similar 

suppression in the C2 case can also come from CiVe- Thus, the actual 6677 rate depends 

strongly also on the mass hierarchy between my^, mg^ and rn^j^. 

Maximizing the Higgsino component of A^i and minimizing that of leads us to the 
conclusion that A''2 is mostly photino and A^i is mostly Higgsino, analogous to the selectron 
interpretation. However, differences do exist between the chargino interpretation and cl models 
in the selectron interpretation. For example, one needs 100 < mg^ ^ 137 GeV in the selectron 
interpretation, while in the chargino interpretation one only needs at least one of e^,, Uf> heavier 
than roughly 60 GeV but lighter than at least one of the charginos. Of course, additional 
constraints on mg^, rriy^ are present, due to the particularly complicated decay chain and the 
large radiative neutralino branching ratio needed. The right selectron enters the 3-body decay 
A''2 — > A^ie''~e~, but if its mass is moderately large then the decay cannot be enhanced. Squark 
masses are relatively unconstrained, although lighter squark masses increase the C2C2 cross 
section, but decrease the radiative branching ratio. 

The absence of 6677 + event kinematical solutions with chargino masses less than 
95 GeV implies that to construct a chargino interpretation that at least possibly satisfies the 
kinematics one should conservatively choose to search only for models with > 95 GeV. 
Restricting to M2 , ^ and tan (3 values roughly in the allowed ranges singled out in the selectron 
interpretation, one finds a rough upper limit of 400, 50 and 1200 fb for the cross section of 
CiCi, C1C2 and C2C2 production respectively. Given at best a useful branching ratio of about 
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5%, then a C1C2 interpretation (alone) can be excluded. For CiCi and C2C2 production, the 
6677 signal could be up to roughly 20 and 60 fb, therefore the lower bound on the radiative 
neutralino decay branching ratio is 50% and 30% respectively, to pass ^ = 5 fb cut used in 
the selectron interpretation. The CiCi cross section drops rapidly as m^^ is increased, and 
it appears not to give a sizeable 6677 + signal when > 110 GeV. Alternatively, the 
C2C2 cross section can still be large, and give a sizeable 6677 signal for ~ 150 GeV (if 
rriq ~ 250 GeV). In practice, this sets rough upper limits for M2 and which determine 
the chargino masses. Further, our analysis of the 6677 + event kinematics in the chargino 
interpretation gives an indication that large mass differences ( > 30 GeV) between Cj and A'^2 
may be required to reconstruct the 6677 + -^^ event. In the i = 1 case this is very difficult, if not 
impossible, given all the other constraints. Thus, we conclude that sizeable 6677 + signals 
can probably only be achieved from (7^(7^ production, with the decay chain C2 i^e^Li^ 6-/V2) 
or C2 6t'e(^ ^€^2), followed by N2 Ni'j. This appears to happen only in a region of the 
parameter space similar to the selectron interpretation. 

A few final remarks on model building are in order. The sneutrino always plays a role 
when the mass hierarchy > me^(> m-Pe) exists, and as a consequence the 6677 signal is 
depleted from C — > t'ee since De tends to have comparable branching ratio into A'^i and A'^2- 
Further, if uii^^ < m^^, then a 2-body decay opens for A''2 — > t'e^^ej which often suppresses the 
radiative decay branching ratio. Also, a sneutrino mass larger than implies a possibly large 
branching ratio for ^ Cie. To ensure sufficient phase space for the decay C2 — > e^fg and 
to have the masses fall in regions where we found kinematical solutions, the mass difference 

— mg^ ^ 0(10) GeV probably should be enforced. The selectron also must be larger 
than ^TT-^^ least ~ 20 GeV for analogous reasons, but not larger than m^^ otherwise the 

branching ratio for cl will be dominated by — > Ciz^e- It is clear that maintaining such a mass 
hierarchy between rn^^, rrigj^, rrii^^, rrij^^, m^_^, is considerably more difficult than in the 

selectron interpretation, and to some extent a fine-tuning of the masses of the particles involved 
is always required. Also, the relevant branching ratio is always small and never exceeds a few 
percent while in the selectron interpretation it can in principle reach 100%. All of these facts 
seem to render a chargino interpretation problematic (in stark contrast to a scenario with the 
gravitino as the LSP Q). 

A. 3 Chargino interpretation — an example 

We searched our model samples compatible with a selectron interpretation of the 6677+^^7^ 
event for cases where C2C2 production could yield an additional 6677 signal. We found several 
tens of candidate models: some in the cr samples, and a few in the ei sample. However, the 
general kinematical requirements for a chargino interpretation of the 6677 + event slightly 
favor the models, which are located roughly in Region 2 (according to the classification of 
Sec. ^^). Such models could give rise to a 6677 signal with the kinematical characteristics of 
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the event, from simultaneously and C2 pair production, although the C2 signal is generally 
below 6 fb. We report one model as an example of the above: Mi = 65 GeV, M2 — Mz, 
^ = -53 GeV tan/? = 2, mg^ = 110 GeV, mg^ = 350 GeV, m£> = 90 GeV, m^^ = 150 GeV, 
« rriq = 250 GeV. The neutralino masses f^fj^^^^ = 65, 70, 96, 137 GeV, and the chargino 
masses rn^^ ^ = 72, 137 GeV. The C2C2 production cross section at the Tevatron is 380 fb, 
while the 6^6^ cross section is 13 fb. The B{N2 Nij) = 81%, the B{C2 t'ec) = 17%, the 
B{C2 eive) = 16%, the B{eL ^ N2e) ~ 100%, and B{ve ^ Cie) = 77%. The 6677 rate is 
roughly 6 fb from only chargino production, and so is slightly above the ^ = 5 fb cut imposed 
in the selectron interpretation. It is worthwhile to remark on how sensitive the 6677 rate is to 
a change in the masses. For example, one can attempt to raise the 6677 rate from chargino 
production by slightly reducing the ei mass in such a way to get a sneutrino lighter than the 
Ci, and gain the additional signal from C2 decays into on-shell sneutrinos and sneutrino decays 
into A''2. This would require mg^ ^ 96 GeV, although the modified model would appear to 
be farther from the region of masses satisfying the 6677 + event kinematics. However, 
the radiative neutralino decay branching ratio drops quite sensitively when the already light 
slepton masses are further reduced. Thus, constructing models in the chargino interpretation is 
somewhat difficult, and it is not obvious how one ought to perturb around any given model to 
increase the 6677 rate. However, we did find some models with interesting characteristics, as 
shown above. A more in-depth analysis is necessary to determine if the chargino interpretation 
is tenable, and if so the ranges of the parameters needed. 



Appendix B: Sample Models 



Here four sample models from the set used in the selectron interpretation are provided 
in Tables 12 and Input parameters and calculated masses are given, along with many 
branching ratios and cross sections. Notice that the four models' input parameters are similar 
(except for the slepton and stop masses) , but the cross sections for both the 6677 + event 
and associated phenomenology are quite different. 
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Model parameters 



bl model 



e/f mo 



del 



Ml, M2 
/X, tan/3 
rriA, rriq = mi 
mi,, Oi 



64.7 , 64.3 
-37.0 , 1.18 

200 , 500 
204 , -0.342 
105 , 272 , 99.6 



74.4 , 77.6 
-38.3 , 1.11 

400 , 500 
487 , -0.123 
391 , 104 , 390 



(iVi|iJ(,)2,(iV2|7)2 
mh,mH,mH±,oih 

B{N2 7Vi7) 

B{N3 ^ l+l-),B{N3 ^ i^u),B{N3_-^ qq) 

B{Na ^uV + w), B{Ni Ij + hi) 

B{Ni^lRl + lRl) 

B{Ci^ mv),B{Ci^ Niqq') 

B{C2^vI),B{C2^1li^) 

B{C2 ^ Nilv),B{C2 ^ Niqq') 

B{eL^ N2e),B{eL-^Cu,) 

B{eR^ N2e),B{eR^ Nic) 

B{i)e^ N3U^),B{i>e^Ce) 



79.6 , 110 
36.6 , 64.6 , 90.5 , 118 
0.997 , 1.000 
70.2 , 229 , 216 , -0.825 
13.2 
0.98 

0.10 , 0.22 , 0.67 
0.83 , 0.13 



0.34 
0.66 
0.02 
0.91 
0.81 
0.08 



0.66 
0.28 
0.03 

0.07 
0.14 
0.90 



78.9 , 119 
38.2 , 75.1 , 88.5 , 127 
0.999 , 0.999 
69.2 , 415 , 408 , -0.765 
6.6 
0.94 

0.10 , 0.20 , 0.69 

0.80 
0.33 , 0.67 

0.33 , 0.66 
0.30 , 0.59 

0.98 , - 
0.10 , 0.61 



LEP161 cross sections: 
a(iViiV3),(T(CiCi) 

a{N2N2),<j{N2N3) 

inclusive a{2l + X), (7(77 + X) 



2010 , 405 
191 , 123 
276 , 184 



2130 , 1320 
40,- 
365 , 36 



LEP190 cross sections: 

a{NiN3),cTiNiNA) 

a{N2N2),'j{N2N3) 

a{CiCi),a{CiC2) 

inclusive a{2l + X), a{'-f^ + X) 

inclusive a{l'y + X), a{llj + X) 



1450 , 89 
342 , 243 

1080 , 167 
473 , 331 
115 , 73 



1530 , 49 
199 , 164 

2760 , - 
529 , 177 
60 , 59 



Tevatron cross sections: 
o-(eLeL),cr(ei{eK) 
a{i>ei>e),(^{eLi>e) 
a{NiN3),cr{CiCi),a{C2C2) 
a{CiNi),a{CiN2),cT{CiN3) 
cT{C2N2),a{C2N3),a{C2N4) 
inclusive a(2l + X), (7(77 + X) 
inclusive a{lj + X), a{2l^ + X) 
inclusive q-(Z77 + X), a{3l + X) 



16.5 , - 
18.5 , 45.0 

1180 , 907 , 552 
2690 , 113 , 840 
324 , 28 , 332 
1700 , 174 
954 , 714 
171 , 892 



-,7.9 

1270 , 887 , 415 
2710 , 55 , 915 
190 , 8.4 , 241 
631 , 24 
318 , 237 
22 , 101 



Table 12: Two sample models in the selection interpretation. All masses are in GeV, all cr( 
sections are in fb. Only the largest branching ratios and cross sections are displayed. I 
summed over e, fi, and r in the branching ratios and inclusive cross sections (which have 
detector efficiencies included). In the branching ratios C refers to a sum over Ci and C2- 
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Model parameters 


+ Gfl, model 


cr model (with light ti) 


Ml, M2 


70.2 , 76.2 


76.5 , 77.0 


/i, tan/3 


-48.8 , 1.26 


-38.9 , 1.39 


ruA, rriq = uif^ 


200 , 500 


400 , 2000 


Til ~ H~ 


488 n 263 


50 , 7r/2 


nil , mr , rno 


119 , 121 , 113 


439 , 105 , 437 




84.8 , 118 


75.2 , 121 


iVi ' iV2 " iV3 ' Jv4 


47.8 , 71.5 , 96.8 , 124 


37.4 , 76.6 , 88.4 , 128 


(iVi|fffa)2,(^2l7)^ 

\ J- 1 1 \ ^\ 1 1 


0.990 , 0.998 


0.988 , 0.999 


l^hi 7 


67.8 , 227 , 216 , -0.792 


59.1 , 411 , 408 , -0.651 




10.2 


5.1 


S(7V2 ^ iVi7) 


0.92 


0.86 


6(7V3 ^ l + l-),B{N3 W),6(7V3 ^ gg) 


0.10 , 0.22 , 0.67 


0.10 , 0.20 , 0.68 


B( Na^ vV + lyv) B(Na ^ItI + It I) 


0.85 , 0.05 




t3( AT ^ T 1 ^ 1 1\ 


U.Ui 


U. / 4 


D(Gi — > iVi/i/j,D(Oi lyiqq ),o[Ci ~* lib) 


U.o4 , U.OD , - 


U.UU , - , i.UU 


B[Ci vl),B{Ci tib) 


0.78 , - 


- , 0.98 


0(62 — > Nilv),B{C2 —> Niqq ) 


O.Oo , 0.11 


0.01 , 0.01 


YJ ( ~ 1\T \ YD 1 " \ 

B{eL -> N2e),B[eL Cve) 


0.94 , 0.03 


0.30 , 0.59 


B{eR N2e),B{eR N^e) 


0.97 , - 


0.96 , - 


B{i)e N3Ve),B{Ue ^ Ce) 


0.10 , 0.86 


0.10 , 0.62 


LEP161 cross sections: 






a{NiN^),a{CiCi) 


1500 , - 


2100 , 2680 




120 , - , - 


23 , - , 850 


inclusive a{2l + X), (7(77 + X) 


157 , 100 


215 , 17 


LEP190 cross sections: 






a(7ViiV3),a(iVi7V4) 


1360 , 24 


1500 , 41 


a{N2N2),aiN2N3) 


355 , 227 


169 , 150 


aiCiCi),aiCiC2),aihil) 


880 , ~ , - 


3110 , " , 760 


inclusive a{2l + X), (7(77 + X) 


302 , 299 


254 , 125 


inclusive cr(/7 + X), ayll^ + X) 


56 , 51 


78 , 78 


Tevatron cross sections: 






cr(eLeL),o-(eHefl) 


9.4 , 4.0 


-,7.5 




10.5 , 24.6 




a{NiN:i),a{CiCi),a{C2C2) 


688 , 681 , 434 


1270 , 1140 , 298 


a{CiNi),a{CiN2),a{CiN3) 


1590 , 86 , 575 


3430 , 128 , 974 


<j{C2N2),a{C2N3),a{C2N4) 


189 , 29 , 259 


218 , 43 , 283 


inclusive a{2l + X), 17(77 + ^) 


1190 , 43 


178 , 16 


inclusive allj + X), a{2l-f + X) 


369 , 279 


50 , 48 


inclusive (7(^77 + X), a{3l + X) 


39 , 654 


16 , 7.5 



Table 13: As in Fig. 12, but for an cl + e_R model, and a model with a light stop. Note that 
a X smxLs over both cl and cr contributions for the eL + e-R model. 
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